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In  the  preceding  Progress  Report,  there  was  considerable  discussion  of  the 
theory  of  energy  bands;  we  spoke  of  the  use  of  a tight  binding  approximation  as  an  in- 
ter polation  method  for  handling  the  problem  of  energy  bands  at  an  arbitrary  point  in 
the  Brillouin  zone,  and  the  orthogonalized  plane  wave  method  for  getting  accurate  re- 
sults at  certain  points  of  special  symmetry  in  the  zone.  Our  thinking  along  these  pnes 
has  progressed  a good  deal  in  the  intervening  time-  In  the  first  place,  in  prenaratiori 
for  the  course  which  I am  proposing  to  give  during  the  fall  on  the  energy  levels  of  im- 
purities in  crystals,  I have  been  thinking  further  about  the  nature  of  Wannier  functions, 
and  about  the  relation  of  the  simplified  tight  binding  approximation  to  such  problems. 
Dr.  Koster  has  been  assisting,  and  Dr.  Schweinler  has  been  making  valuable  sugges- 
tions in  connection  with  these  questions,  and  a short  summary  of  some  of  our  results 
is  included  in  this  Report,  though  they  will  be  written  up  m more  detail  in  the  Technical 
Report  No.  5,  which  will  include  the  material  to  be  given  during  this  course.  Some  of 
the  points  which  we  are  finding  in  this  study  begin  to  have  a rather  close  relationship 
to  the  actual  behavior  of  impurity  levels  in  germanium  and  silicon,  such  as  the  possi- 
bility of  two  quite  different  sorts  of  impurity  levels,  hydrogen-like  ones  of  the  sort 
often  discussed,  and  also  much  deeper-lying  levels  with  the  wave  functions  more  con- 
centrated around  the  impurity  atoms.  In  this  discussion  of  impurity  levels,  we  had 
the  benefit  of  a visit  from  Professor  H.  Jones  of  the  Imperial  College  of  Science  and 
Technology,  Dr,  J.  Friedel  of  the  Ecole  des  Mines,  and  Professor  C.  Kittel  of  the 
University  of  California,  during  which  the  current  state  of  such  theory  was  discussed; 
much  of  the  progress  reported  here,  however,  has  been  made  since  that  visit  during 
July. 

The  other  aspect  of  the  theory  of  energy  bands  in  which  we  have  made  progress 
is  an  improvement  in  the  orthogonalized  plane  wave  method.  Dr.  Parmenter  in  the 
last  Progress  Report  described  a method  which  he  has  been  applying  to  chromium;  in 
thinking  about  that  method,  smd  about  the  method  as  proposed  by  Herring  and  applied 
to  diamond  by  Herman,  it  occurred  to  me  that  there  was  a variation  of  the  existing 
methods  which  might  have  some  advantages  over  any  of  them.  This  consists  in  a 
scheme  to  satisfy  Schrodinger 's  equation  as  accurately  as  possible  within  the  atoms, 
where  we  assume  spherical  symmetry,  amd  to  fit  such  a solution  onto  plane  waves  out- 
side the  atoms,  then  making  a linear  combination  of  such  augmented  plane  waves  to 
get  a final  wave  function.  I had  tried  such  a scheme,  in  1937,  but  the  method  used  then 
was  not  conveniently  adapted  to  calculation,  and  it  has  had  very  little  use.  The  new 
scheme,  however,  avoids  the  difficulties  of  the  old  one,  and  we  hope  that  it  may  prove 
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useful,  and  rapidly  convergent,  '.  he  mathematicc.1  dcvclcpment  as  1 carried  it  through, 
involved  some  rather  intractable  infinite  sums.  However,  when  Mr.  Saffren,  who  has 
newly  joined  the  Group,  looked  into  the  mathematical  formalism,  he  v/as  able  by  a very 
ingenious  application  of  Green's  functions  to  show  how  these  series  could  be  summed, 
leading  finally  to  a very  simple  and  usable  formalism.  We  hope  that  the  method  can 
shortly  be  applied  to  several  crystals.  We  are  fortunate  in  that  Dr.  D.  J.  Howarth,  of 
the  Imperial  College,  has  just  joined  the  Group  for  a visit  of  a year.  He  has  been  work- 
ing with  Jones  on  very  accurate  applications  of  the  cellular  method  to  sodium,  copper, 
and  other  crystals.  We  are  hoping  that  in  collaboration  with  him  we  can  work  out  a pro- 
gram of  applying  the  new  method  to  these  same  crystals,  so  as  to  compare  the  two 
rather  different  approaches.  The  results  of  accurate  calculations  of  energy  bands  are 
so  complicated  that  there  is  a very  natural  question  as  to  how  close  the  final  res  iUs 
are  to  real  solutions  of  the  problem.  Our  confidence  would  oc  greatly  improved  by  get- 
ting results  by  two  methods  which  checked  in  their  details.  At  the  same  time,  such  a 
comparison  should  show  which  method  is  easier  to  apply. 

We  come  next  to  the  many  molecular  calculations,  principally  involving  con- 
figuration interaction,  which  are  going  on.  In  this  connections,  we  include  a report 
from  Dr.  R.  McWeeny,  of  King's  College,  University  of  Durham,  located  at  Newcastle- 
upon-Tyne,  who  has  just  joined  the  Group  for  a year's  visit.  Dr.  McWeeny  has  been 
investigating  problems  of  configuration  interaction,  and  comparing  the  results  with 
those  obtained  from  simple  wave  functions  constructed  from  molecular  orbitals,  and 
also  with  Heitler -London  wave  functions,  constructed  according  to  the  standard  Heitler- 
London  method,  but  with  the  calculations  carried  out  properly,  involving  all  the  terms 
usua’ly  omitted.  His  actual  calculations  have  been  made  using  orthogonalized  orbitras, 
which  we  believe  is  the  most  practicable  method  for  complicated  systems.  The  inte:*- 
esting  result  which  he  finds  is  that  the  Heitler-London  wave  functions  are  rather  close 
to  the  best  wave  functions  derived  from  configuration  interaction,  much  closer  to  them 
than  are  the  molecular  orbital  wave  functions.  This  is  what  wc  should  suspect  from 
the  example  of  hydrogen,  where  it  is  certainly  true.  The  importance  of  studies  like 
that  of  McWeeny  is  that  if  we  find  that  a wave  function  like  the  Heitler-London  one  really 
is  rather  accurate,  we  can  then  express  this  wave  function  as  a combination  of  a num- 
ber of  configurations  formed  from  orthogonalized  orbitals  (so  as  to  put  the  calculations 
in  a form  which  can  readily  be  carried  out),  and  it  is  likely  to  be  a rather  good  first 
guess  at  the  combination  of  configurations  which  will  minimize  the  energy. 

Any  such  guide  to  the  correct  configuration  interaction  will  be  of  great  value  in 
more  complicated  problems,  whereas  we  have  seen  in  earlier  reports  a complete  con- 
figuration interaction  is  too  complicated  to  undertake.  Such  a use  of  the  Heitler-London 
method  is  cf  course  completely  different  from  the  misuse  of  the  method  which  has  so 
often  been  made,  in  which  many  integrals  are  neglected  in  the  interest  of  simplification. 
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Such  approximations  arc  not  involved  in  the  present  method,  which  gets  at  the  Heitler- 
I ondon  results  through  the  use  of  configuration  interaction  methods  with  orthogonalized 
orbitals.  These  interesting  results  of  McWeeny  are  included  in  this  Progress  Report, 
though  they  were  obtained  before  he  joined  the  Group,  for  two  reasons;  because  they 
contribute  so  directly  to  the  line  of  thought  on  configuration  interaction  which  we  have 
been  following,  and  because  he  will  continue  to  work  along  the  same  line,  and  make 
contributions  to  future  progress  reports  based  on  the  material  presented  here. 

The  calculation  on  antiferromagnetism  which  Dr.  Pratt  reports  has  a good  deal 
of  resemblance  to  that  of  McWeeny,  in  a way.  In  the  previous  Progress  Report  he  had 
indicated  the  difficulty  of  getting  a model  of  a molecule  resembling  (MnOMn)^^  which 
did  not  lead  to  ferromagnetism  rather  than  antiferromagnetism,  and  we  suspected  that 
a configuration  interaction  was  required  to  lead  to  the  observed  results.  Since  he  did 
not  seem  to  be  able  to  choose  the  correct  configurations  to  lead  to  the  desired  results, 
when  expressed  in  terms  of  orthogonalized  orbitals,  I suggested  that  he  set  up  the  prob- 
lem in  terms  of  non-orthogonalized  orbitals,  following  a Heitler-London  model,  using 
a simplified  four -electron  case  which  might  lead  to  the  antiferromagnetic  sort  of  action, 
and  carrying  the  configuration  interaction  through  completely  with  the  non-orthogo- 
nalized orbitals.  This  use  of  non-orthogonalized  orbitals  in  a Heitler-London  calcula- 
tion, as  I have  pointed  out  before,  very  rapidly  becomes  entirely  impossible  as  the 
number  of  electrons  increases,  but  for  four  electrons  it  is  still  possible,  though  diffi- 
cult. Fortunately  this  method  seems  to  have  worked;  and  now,  by  expanding  the  final 
wave  function  in  terms  of  orthogonalized  orbitals,  Pratt  has  found  which  configurations 
are  required  to  get  the  desired  results,  and  which  he  had  failed  to  include  in  his  ear- 
lier work.  It  is  interesting  that  this  exact  Heitler-London  calculation  can  be  compared 
with  the  simplified  type  of  approximate  use  of  the  method  which  we  feel  is  unreliable. 

In  this  case,  even  though  the  overlap  integral  is  only  of  the  order  of  0.  04,  an  error  of 
about  fifty  pcrcern  is  made  by  using  the  usual  Dirac  vector -model  methods,  neglecting 
o\'erlap  integrals  and  all  exchange  integrals  except  the  two-center  integrals.  This 
ce.rtainly  should  put  one  on  his  guard  against  using  that  method  for  any  work  in  which 
it  is  regarded  as  an  accurate  method  of  calculation. 

Pratt's  results  are  interesting  physically,  as  well  as  pointing  the  way  toward 
better  methods  of  calculation.  They  show,  more  or  less  as  in  the  discussion  of  the 
preceding  Progress  Report,  that  the  physically  important  effect  is  a distortion  of  the 
oxygen  atom  (here  replaced  for  simplification  by  a helium  atom)  by  the  two  adjacent 
magnetic  atoms.  The  detailed  mechanism,  however,  is  not  that  described  in  the  pre- 
ceding Pr.igress  Report,  but  instead  is  the  process  described  as  super -exchange,  by 
Anderson  and  Van  Vleck;  we  feel,  in  other  words,  that  their  mechanism  is  correct, 
even  though  we  do  not  feel  that  the  method  used  by  them  to  compute  it  is  reliable. 

There  are  two  possible  types  of  distortion,  or  polarization,  of  the  oxygen  atom,  either 
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of  which  could  lead  to  a stabilization  of  the  singlet  state  compared  to  the  triplet.  In 
either  mechan’sm,  the  electrons  of  plus  spin  in  the  oxygen  will  be  displaced  toward 
one  neighboring  atom,  those  of  minus  spin  toward  the  other,  and  this  resulting  defor- 
mation can  lower  the  energy,  being  a form  of  polarization.  The  mechanism  described 
in  the  preceding  Progress  Report  was  one  in  which  the  electron  of  plus  spin  from  the 
oxygen  wishes  to  be  attracted  toward  the  ion  with  plus  spin,  on  account  of  Hund’s  rule. 
The  opposite  mechanism,  the  superexchange  mechanism,  is  that  in  which  it  is  attracted 
to  the  ion  with  minus  spin,  on  account  of  an  incipient  covalent  binding.  One  might  well 
expect  that  the  first  would  be  the  likely  mechanism  in  case  the  magnetic  ion  had  an 
only  partially  filled  d shell,  so  that  the  oxygen  electron  attracted  to  it  could  enter 
this  shell  with  the  same  spin  without  contradicting  the  exclusion  principle;  whereas 
the  second  mechanism  would  apply  when  all  the  orbitals  of  one  spin  in  the  d shell  were 
filled,  all  those  of  the  other  spin  empty,  as  in  MnO,  and  as  in  the  He^  model  which 
Pratt  used.  In  either  case,  Pratt's  present  contribution  points  the  way  to  a quantita- 
tive method  of  handling  the  effect. 

Another  magnetic  calculation,  pointing  out  the  complications  of  the  problem 
and  the  dangers  of  trusting  to  over-simplified  exchange  integrals,  is  Kikuchi's  calcula- 
tion of  the  interaction  of  two  p electrons,  to  see  if  they  can  have  a ferromagnetic  in- 
teraction, as  in  This  very  simple  problem,  as  Kikuchi  shows,  leads  to  a four-by- 
four  secular  equation,  and  no  one  who  has  looked  at  his  figure,  illustrating  the  way  in 
which  the  energy  levels  of  singlet  and  triplet  states  vary  with  internuclear  distance  in 
this  very  simple  problem,  can  doubt  the  extreme  complexity  of  configuration  interac- 
tions, and  the  dangers  of  oversimplifying  them,  it  is  very  interesting  that  Kikuchi 
found  a magnetic  ground  state  over  a range  of  internuclear  distances,  though  the  singlet 
lies  lower  at  smaller  internuclear  distances,  as  it  should  from  general  grounds.  But 
the  behavior  at  large  distances  is  so  complicated  that  it  is  clear  that  no  method  based 
on  simple  exchange  integrals,  amounting  to  an  extrapolation  inward  from  infinite  sepa- 
ration, could  be  reliable  in  such  a case.  Kikuchi's  case  is  far  enough  from  actual 
magnetic  materials  so  that  it  can  hardly  be  applied  to  an  actual  case,  but  it  is  very 
valuable  to  see  how  these  problems  work  out  in  simple  cases.  Dr.  Kikuchi  has  left 
the  Group  since  the  last  Progress  Report,  to  take  a position  at  the  University  of  Chicago, 
and  his  absence  will  be  seriously  felt. 

Still  another  contribution  to  the  study  of  configuration  interaction  is  the  result 
which  Mr.  Callen  has  obtained  in  his  study  of  configuration  interaction  in  hydrogen. 
Coulson  a number  of  years  ago  showed  that  by  getting  a really  good  symmetric  molecu- 
lar orbital  for  the  ground  state  of  hydrogen,  the  molecular  orbital  energy  was  im- 
proved a good  deal  as  compared  to  the  calculations  made  using  the  LCAO  MO's.  Callen 
was  interested  in  setting  up  an  equally  good  antisymmetric  molecular  orbital,  and 
then  finding  the  result  of  a configuration  interaction  between  just  the  two  configurations 
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in  which  both  electrons  are  in  the  symmetric  state,  or  both  in  the  antisymmetric  state. 
It  is  this  type  of  configuration  interaction  using  LCAO  MO's  which  was  discussed  at 
length  in  I'echnical  Report  No.  3,  and  which  includes  the  Heitler-London  calculation 
as  a special  case,  but  which  does  not  lead  to  very  good  numerical  results.  Callen  has 
carried  through  this  calculation,  and  his  present  result  for  energy  lies  among  the 
better  approximations  to  the  hydrogen  problem,  aside  of  course  from  the  James  and 
Coolidge  result,  which  is  of  quite  different  order  of  accuracy.  He  expects  to  be  able 
to  improve  his  energy  somewhat  by  further  varying  his  antisymmetric  wave  function. 

The  various  other  pieces  of  work  under  way  in  the  Group  are  making  progress, 
as  described  in  the  Report,  though  they  have  not  reached  a point  of  interesting  results. 
There  have  been  certaii.  changes  in  the  personnel  of  the  Group,  since  the  last  Report, 
in  addition  to  those  already  noted:  Dr.  Barrett  has  completed  the  term  of  his  fellow- 
ship, and  has  left  to  take  a position  at  the  North  Carolina  State  College,  and  Mr.  Hunt 
has  entered  the  Army. 


1.  WAVE  FUNCTIONS  FOR  IMPURITY  LEVELS 


It  is  well  known  that  the  best  method  for  handling  the  wave  functions  of  per- 
turbed periodic  lattices  is  by  means  of  the  W'annier  functions.  We  can  treat  Schro- 
dinger's  equation  in  such  a case  by  making  a superposition  of  Wannier  functions  on 

different  atoms,  and  in  different  energy  bands.  If  a (r  - R.)  is  a Wannier  function  in 

tp  ^ ^ 

the  n band,  located  on  the  atom  in  the  unit  cell  at  vector  position  R.,  then  such  a 

superposition  is  ^(n,  R.)  U^(H.)  Schrodinger ’s  equation  in  the  presence 

of  a perturbative  potential  V,  superposed  on  the  periodic  potential  for  which  the  ori- 
ginal energy  bands  are  calculated,  reduces  to  a set  of  linear  equations  for  the  U'.s, 
which  can  be  written  in  the  form 


‘(m,  R)H  (R.,R.)U  (R  ) = EU  (R.) 
j'  nm'  i'  j'  m j n r 


(1.1) 


The  matrix  components  of  the  Hamiltonian,  ^j)  composed  of  the  contribu- 
tions from  the  periodic  potential,  which  can  be  written  ^ (R.  - R.)  6 , and  the  ma- 

n 1 J nm 

trix  components  perturbation.  We  see  that  th . periodic  potential 

has  no  matrix  components  between  Wannier  functions  in  different  bands,  and  that  the 
matrix  components  between  Wannier  functions  on  different  atoms  but  in  the  same  band 
depend  only  on  the  displacement  from  one  atom  to  the  other. 

The  reason  why  the  Wannier  function  method  has  appealed  so  much  to  the  physi- 
cal intuition  is  that  the  function  U^(R^)  has  itself  many  of  the  aspects  of  a solution  of 
Schrodinger 's  equation.  We  may  replace  this  set  of  coefficients,  defined  for  discrete 
R.'s,  by  a continuous  function  U^(r).  We  may  then  approximate  the  secular  equation 
for  the  U^^(R.)'s,  which  is  a set  of  difference  equations,  by  differential  equations,  and 
if  the  U's  vary  slowly  enough  from  atom  to  atom,  these  differential  equations  reduce 
to  second  order  differential  equations,  which  are  equivalent  to  Schrodinger  equations 
for  a particle  with  an  effective  mass  moving  in  the  perturbative  field.  These  approxi- 
mations are  not  always  very  valid,  however,  ahd  we  have  been  looking  into  the  straight- 
forward problem  of  solving  the  difference  equations,  which  is  a much  more  valid  pro- 
cedure. 


As  an  example  of  a case  where  the  differential  equation  is  definitely  inappli- 
cable, but  where  the  difference  equation  methexi  can  be  carried  through,  we  may  con- 
.sider  the  problem  of  a perturbing  atom  in  a three-dimensional  crystal,  in  a case  where 
we  can  get  a good  approximation  by  using  the  Wannier  function  of  only  one  band.  In 
particular,  w'c  have  considered  that  ^j)  reduced  to  a single  term,  for  n = m 

(since  we  consider  only  one  band),  and  for  = R^  = 0.  In  other  words,  we  assume 
that  the  perturbing  potential  is  so  concentrated  that  it  is  felt  only  on  the  central  atom. 
This  is  not  a bad  approximation  for  the  effect  of  an  impurity  atom  in  a metal,  where 
the  perturbative  potential  is  known  to  be  screened  very  strongly  from  even  close 
neighbors.  We  have  further  simplified  the  problem  to  the  case  of  a simple  cubic  lattice. 
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with  matrix  components  ■ ^j)  connected  with  only  nearest  neighbor  interactions. 

In  such  a case  the  difference  equation  can  be  solved  exactly,  by  an  adaptation  of  a 
method  used  by  McCrae  and  Whipple,  which  waS  pointed  out  to  us  by  Dr.  Schweinler. 
We  shall  not  reproduce  the  solution  here;  it  is  somev/hat  complicated,  and  will  be 
given  in  detail  in  Technical  Report  No.  5.  We  can  state  the  general  results,  however. 

We  can  get  a guide  as  to  what  to  expect  from  the  corresponding  free  electron 
problem,  representing  the  differential  equation  which  approximates  the  difference 
equation.  This  is  the  problem  of  a free  electron  with  a deep  spherical  potential  well 
near  the  origin,  but  otherwise  in  a constant  potential.  We  know  that  this  three-dimen- 
sional potential  well  problem,  unlike  the  corresponding  one-dimensional  problem,  does 
not  liave  any  solutions  of  the  form  of  discrete  energy  levels  unless  the  well  exceeds 
certain  critical  dimensions.  The  reason  is  simple.  The  solutions  of  the  three-dimen- 
sional problem  inside  the  potential  well  is  of  the  form  sin  kr/r,  and  outside  of  the  form 
e r.  The  functions  sin  kr  and  e must  join  continuously  and  with  continuous 
slope  at  the  boundary  of  the  well.  This  is  clearly  impossible  unless  the  function  sin 
kr  has  gone  through  at  least  a quarter  wave  length  inside  the  well,  so  that  it  has  risen 
above  its  maximum  and  is  sloping  down  when  it  must  join  the  decreasing  exponential. 

If  the  well  is  just  deep  enough  to  join  onto  an  exponential  with  negligible  y,  the  wave 
function  outside  falls  off  like  l/r,  and  if  it  is  deeper,  the  wave  function  e*  falls 
off  even  faster  than  l/r. 

Our  study  of  the  difference  equation  shows  that  here  too  there  is  a critical  value 
of  the  pertu*'bation  of  the  potential  on  the  central  atom,  which  must  be  exceeded  before 
a discrete  level  appears.  Numerically,  the  perturbative  term  V^  must  be  equal  to 
approximately  a third  of  the  band  width  before  there  is  such  a discrete  state.  By  the 
time  V^  becomes  equal  to  something  more  than  half  the  band  width,  the  wave  function 
of  the  discrete  level  is  falling  off  so  rapidly  that  almost  the  whole  of  the  wave  function 
is  confined  to  the  central  atom.  In  any  case,  so  long  as  there  is  a discrete  state  at 
all,  the  wave  function  falls  off  rapidly  with  r (as  it  falls  off  at  least  as  fast  as  l/r  in 
the  case  of  the  differential  equation),  and  this  is  such  a rapid  variation  with  position 
that  it  is  not  legitimate  in  any  case  to  replace  the  difference  equation  by  a differential 
equation,  for  quantitative  work.  The  differential  equation  is  useful  only  as  a general 
guide  to  the  nature  of  the  results. 

Arguments  of  the  type  we  have  just  been  sketching,  though  worked  out  on  the 
basis  of  the  differential  equation,  have  been  used  by  Friedel^^^  to  discuss  impurity 
levels  in  alloys;  he  makes  much  use  of  the  fact  that  in  some  cases  there  are  discrete 
levels,  and  in  other  cases  not.  Clearly  this  distinction  can  be  important  in  the  con- 
sideration of  the  effect  of  impurity  atoms  on  energy  bands.  It  appears  in  addition, 
however,  that  similar  questions  enter  in  the  problem  of  impurity  levels  in  semicon- 
ductors. The  problem  <>r  perturbed  levels  in  the  presence  of  a Coulomb  potential  is 
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quite  different.  fr'>m  that  of  the  square  well  potential.  We  know  that  the  differential 
equation  has  an  infinite  number  of  discrete  solutions,  rather  than  a finite  number  of 
zero  as  with  the  square  well.  Furthermore,  the  Is  wave  function  falls  off  like  e 
rather  than  as  e r as  wiih  the  square  well.  For  the  problem  of  impurity  levels  in 
a dielectric,  like  germanium,  with  high  dielectric  constants,  it  is  a familiar  fact  that 
the  constant  y is  so  small  that  the  wave  function  falls  off  only  slightly  from  one  atom 
to  its  neighbor.  Thus  in  this  case  it  is  legitimate  to  replace  the  difference  equation  by 
the  differential  equation  in  the  standard  way,  and  that  method  should  give  good  approxi- 
mations. 

There  is  one  exception  to  this  statement,  however.  In  the  Coulomb  case,  we 
assum"!  that  the  significant  matrix  components  of  perturbative  potential,  ^j)> 

come  when  n = m,  i?.  = R^,  and  we  assume  that  the  components  ^^^n^^i’  ^i^  represent 
Coulomb  energies,  as  modified  by  the  dielectric  constant,  on  atoms  at  a distance 
from  the  perturbing  atom.  On  account  of  the  large  dielectric  constant,  these  terms 
are  small.  On  the  other  hand,  the  dielectric  constant  does  not  enter  into  the  perturba- 
tive atom  on  the  central  atom  itself.  The  term  0)  can  be  quite  large.  It  is  of 

course  not  infinite,  as  a Coulomb  potential  would  be  ai  the  origin;  for  it  represents  the 
average  of  a perturbative  potential  of  Coulomb  form  (without  the  dielectric  correction) 
over  the  Wannier  function  on  the  central  atom.  Nevertheless  it  is  a large  term,  which 
is  not  considered  at  all  in  the  usual  formulation  of  the  hydrogen-like  impurity  levels  in 
a semiconductor. 

It  appears  likely  that  this  term  V^(0,  0)  will  act  much  like  the  corresponding 
term  in  the  problem  of  a perturbation  on  a single  atom,  which  we  have  already  men- 
tioned. That  is,  unless  it  exceeds  a critical  value,  it  is  not  likely  to  make  much  modi- 
fication of  the  hydrogenic  wave  function  and  energy  level,  as  usually  considered.  If  it 
exceeds  a critical  value,  however,  it  can  modify  the  character  of  the  wave  functions 
completely,  causing  them  to  be  concentrated  much  more  closely  around  the  perturbing 
atom,  and  can  push  the  discrete  energy  level  far  away  from  tlic  valence  or  conduction 
band.  We  must  look  in  some  cases,  in  other  words,  for  energ;/  levels  midway  between 
valence  and  conduction  bands;  this  may  well  explain  the  origin  of  the  deep  trapping 

(4) 

levels  which  are  often  observed  in  semiconductors. 

In  the  ordinary  case  of  impurity  levels,  it  is  not  legitimate  to  disregard  the 

interaction  of  several  energy  bands,  with  their  effect  on  the  symmetry  of  the  Wannier 

(5) 

functions.  It  has  already  been  pointed  out  ' that  the  Wannier  functions  with  overlapping 
bands  have  very  complicated  symmetry  properties.  In  such  a case,  however,  we  be- 
lieve that  the  most  practicable  approach  is  to  use  the  simplified  tight  binding  approxi- 
mations^^ which  has  been  discussed  in  the  preceding  Progress  Report,  and  in  Technical 
Report  No.  4.  As  an  example,  let  us  take  the  diamond  lattice,  as  in  silicon  or  ger- 
manium. There  it  is  possible  to  set  up  good  approximations  to  the  valence  and  conduc- 
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tion  bands  by  using  eight  types  of  atomic  orbitals,  the  2s,  2p  , 2p  , 2p  orbitals  on 

X y z 

each  of  the  two  atomic  f.ites  in  the  unit  cell.  We  can  now,  to  an  equally  good  approxi- 
mation, build  up  the  wave  functions  of  impurity  levels  by  superposing  the  same  types 
of  atomic  orbitals.  That  is,  we  set  up  a superposition  ^ (n,  R.)  U^(R^) 
where  the  4>'s  are  the  atomic  orbitals,  to  represent  the  wave  function,  where  the  sum- 
mation is  over  these  eight  orbitals.  The  difference  equations  for  determining  the  U's 
are  exactly  of  the  form  of  Eq.  (l),  except  that  the  matrix  components  for  the  oeriodic 
potential,  instead  of  being  ^^\Rj  - iT^)  replaced  by  the  matrix  components 

of  energy  between  the  various  orbitals  set  up  in  the  tight -binding  approximation  in 
Ref.  (6).  If  we  choose  these  matrix  components  to  be  identical  to  those  found  in  Ref. 
(6),  chosen  to  give  go<xl  agi  ecineiu  with  ihe  eoi  l eci  calculation  of  the  energy  bands  at 
points  in  the  Brillouin  zone  where  they  have  been  worked  out  by  the  orthogonalized 
plane  wave  method  or  other  accurate  approximation,  then  the  wave  functions  of  im- 
purity levels  set  up  in  this  way  will  be  entirely  consistent  with  the  theory  of  the 
energy  bands. 

This  procedure  can  be  shown  to  be  equivalent  to  the  following  much  more  com- 
plicated procedure:  to  solve  for  the  energy  bands  (eight,  in  the  case  of  the  diamond 
lattice)  by  the  tight -binding  approximation,  to  compute  Wannier  functions  from  each 
one  of  the  resulting  (eight)  energy  bands,  and  to  use  these  Wannier  functions  for  de- 
termining the  impurity  wave  functions,  with  use  of  interactions  between  all  these  Wan- 
nier functions.  By  using  the  atomic  orbitals  directly,  we  can  see  much  more  clearly 
the  nature  of  the  resulting  impurity  wave  function.  Each  wave  function  will  have  a 
definite  symmetry,  transforming  according  to  an  irreducible  representation  of  the 
point  symmetry  group  of  the  crystal.  We  can  then  see  at  once  that  on  each  atom  we 
must  use  combinations  of  the  appropriate  atomic  orbitals  so  that  the  whole  wave  func- 
tion will  have  the  correct  symmetry.  This  can  often  give  us  enough  infoi  mation  so 
that  we  can  predict  in  advance  that  many  of  the  U^(^^)’s  must  be  zero,  or  have  certain 
relations  to  each  other.  As  an  example  of  the  sort  of  thing  to  expect,  we  should  ex- 
pect that  an  impurity  level  of  s-like  symmetry,  constructed  largely  out  of  p atomic 
orbitals,  would  have  the  axes  of  these  p orbitals  pointing  approximately  radially  in- 
ward or  outward  from  the  central  atom.  This  is  similar  to  the  construction  of  a 
Wannier  function  of  s-like  symmetry,  from  p orbitals,  discussed  in  Ref.  5. 
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Z.  AN  AUGMENTED  PLANE  WAVE  METHOD  FOR  THE 

PERIODIC  POTENTIAL  PROBLEM 

The  cellular  method,  and  the  orthogonalized  plane  wave  method,  form  two  ap- 
proximations to  the  solution  of  the  nof-odic  potential  problem  which  are  so  well  known 
that  they  need  no  comment.  Both  of  them  take  advantage  of  the  fact  tlirt  the  poten- 
tial is  approximately  spherically  symmetric  about  each  atom.  In  fact,  we  can  get  a 
good  enough  approximation  to  the  periodic  potential  in  a real  crystal  to  serve  as  a 
starting  point,  by  assuming  that  it  is  spherically  symmetrical  within  atomic  spheres, 
constant  between  the  spheres;  we  may  best,  in  such  a case,  take  the  spheres  as  large 
as  possible,  preferably  touching.  Then  it  is  a natural  hope  that  somehow  we  can  build 
up  good  solution*:  of  thc  composite  SchrSdiriger  piuhlem,  since  we  Know  how  lo  solve 
this  problem  exactly  within  the  atomic  spheres  where  the  potential  is  spherically  sym- 
metric, and  also  in  a region  of  constant  potential  such  as  the  region  outside  the  atoms. 
One  of  the  present  writers^^^  tried  in  1937  to  set  up  a solution  in  the  form  of  a super- 
position of  functions  each  of  which  was  a plane  wave  between  the  spheres,  joined  con- 
tinuously onto  a solution  of  Schrodinger's  equation  within  each  of  the  spheres.  Func- 
tions of  this  type  had  some  resemblances  to  the  orthogonalized  plane  waves  introduced 
later  by  Herring  for  the  same  purpose,  in  that  in  each  case  it  was  recognized  that  the 
difficulty  with  using  ordinary  plane  waves  to  solve  Schrodinger's  equation  was  that  they 
formed  such  very  poor  approximations  near  the  atoms.  The  method  mentioned  in  Ref. 
(2)  never  had  much  application,  because  certain  features  made  it  very  awkward  and 
impracticable  to  use.  The  present  note  outlines  a new  scheme,  with  the  same  aim, 
which  gives  promise  of  being  much  more  practical. 

The  whole  problem  is  to  take  a single  plane  wave  in  the  region  between  the 
spheres,  and  match  this  as  well  as  possible  to  a solution  of  the  spherical  Schrodinger 
probleui  within  each  sphere.  The  matching  cannot  be  perfect;  it  is  well  known  from 
scattering  theory  that  a plane  wave  outside  a scattering  atom  can  be  matched  to  a solu- 
tion inside  only  if  we  include  scattered  waves  outside  the  atom.  Our  aim  rather  is  to 
do  the  best  matching  job  possible,  and  then  take  care  of  the  scattering  problem  by 
superposing  a number  of  such  plane  waves,  which  we  may  call  augmented  plane  waves, 
since  they  contain  an  additional  spherical  contribution  within  each  atom.  Our  aug- 
mented plane  waves  will  be  similar  to  Herring's  orthogonalized  plane  waves,  only  we 
hope  that  by  paying  very  close  attention  to  satisfying  Schrodinger's  equation  within  the 
atoms  we  can  make  them  better  initial  approximations  to  exact  solutions,  and  conse- 
quently can  make  the  process  of  combining  a number  of  such  augmented  plane  waves 
more  rapidly  convergent  than  by  Herring's  method.  Herman  and  Callaway,  in  their 
study  of  diamond  and  germ*inium,  have  shown  that  the  process  converges  sufficiently 
fast  to  be  usable.  The  reason  why  we  prefer  to  handle  the  scattering  problem  by  a 
superposition  of  plane  waves,  rather  than  of  scattered  wavelets,  as  in  the  problem  of 
scattering  by  single  atoms,  is  that  using  the  plane  waves  we  can  immediately  satisfy 
the  periodicity  conditions  in  the  crystal,  whereas  if  we  superpose  scattprea  wavelets. 
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as  Korringa  does,  we  have  a major  problem  to  satisfy  these  periodicity  conditions. 

The  matching  procedure  which  is  now  being  proposed  is  the  following.  We 
start  with  a plane  wave  of  definite  propagation  constant  k,  energy  E^,  in  the  region 
outside  the  atoms,  where  the  constant  potential  energy  is  set  equal  to  zero,  so  that  in 
atomic  units  E^  = k^.  Around  the  boundary  of  an  atom,  we  can  expand  this  plane 
wave  in  a series  of  spherical  harmonics  of  angle,  and  spherical  Bessel  functions  of 
r,  in  the  well  known  way.  Inside  the  sphere,  we  now  set  up  Schrodinger 's  equation 
for  the  assumed  spherical  potential,  for  each  1 value.  We  solve  this  equation  using 
the  energy  as  a parameter,  taking  the  solutions  regular  at  the  origin,  and  find  the 
logarithmic  derivative  of  the  function  at  the  surface  of  the  sphere,  as  a function  of 
the  energy.  We  pick  out  those  particular  energies,  and  eigenfunctions,  for  which  the 
logarithmic  derivative  equals  the  logarithmic  derivative  of  the  corresponding  spheri- 
cal Bessel  function  encountered  in  the  expansion  of  our  particular  plane  wave.  If  we 
make  any  finite  sum  of  such  eigenfunctions,  then,  with  coefficients  so  chosen  that  the 
function  joins  continuously  onto  the  plane  wave  at  the  boundary,  we  shall  find  that  the 
derivative  of  the  function  also  joins  smoothly;  for  the  functions  inside  and  outside 
have  the  same  logarithmic  derivative,  or  ratio  of  slope  to  function,  at  the  joining  sur- 
face. We  can  build  up  a function  inside,  then,  joining  smoothly  onto  the  plane  wave, 
in  a very  large  number  of  ways,  since  we  can  give  arbitrary  coefficients  to  all  of  these 
eigenfunctions,  subject  only  to  the  condition  that  the  function  be  continuous  at  the  sur- 
face of  the  sphere. 

The  next  step  in  the  procedure  is  to  set  up  the  expectation  value  of  the  energy, 
for  the  function  which  equals  the  plane  wave  outside  the  sphere,  the  arbitrary  super- 
position of  our  eigenfunctions  within  the  sphere,  and  to  vary  the  arbitrary  parameters 
to  make  the  energy  stationary.  By  so  doing,  we  may  hope  to  get  the  best  approximate 
solution  of  Schrodinger 's  equation,  consistent  with  equalling  the  plane  wave  outside  the 
sphere,  and  joining  smoothly  onto  it  at  the  surface  Fortunately,  this  variation  pro- 
cedure can  be  carried  through  explicitly,  without  having  to  solve  an  unmanageable 
secular  equation.  We  can  set  up  an  explicit  formula  for  the  resulting  augmented  plane 
waves  and  their  energies.  We  find  that,  for  a given  k value,  we  have  solutions  of  low 
energies,  corresponding  to  »he  bound  states  of  the  atoms.  The  energies  are  almost 
exactly  identical  with  those  of  Ihe  free  atoms,  and  the  wave  functions  are  concentrated 
almost  completely  within  the  atoms,  the  amplitudes  of  the  plane  waves  between  the 
atoms  being  negligible.  For  higher  energies,  on  the  conl.rary,  we  have  solutions 
which  are  much  more  like  the  ordinary  orthogonalized  plane  waves;  they  are  of  com- 
parable magnitudes  inside  and  outside  the  atoms,  joining  smoothly  over  the  surface 
of  the  atom.  These  functions  have  some  very  useful  properties.  Functions  of  the 
same  k value  are  automatically  orthogonal  (so  that  we  have  achieved  the  orthogonali- 
zation  which  Herring  uses)  and  have  no  non-diagonal  matrix  components  of  energy 
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between  them,  since  they  are  derived  from  a variation  method.  We  have  non-vanishing 
non-diagonal  matrix  components  of  energy  between  augmented  waves  of  different  k 
values,  however,  and  these  components  are  easy  to  compute,  and  are  the  quantities 
required  in  setting  up  the  linear  superposition  of  augmented  plane  waves  which  is  to 
be  used  as  an  approximation  to  the  true  wave  function. 

The  method  in  this  form  was  worked  out  by  one  of  us  (JCS),  and  is  appearing 
in  a forthcoming  number  of  the  Physical  Review.  It  has  the  disadvantage  that,  though 
the  formulation  is  simple  and  explicit,  still  it  contains  a number  of  infinite  sums,  in 
the  form  of  sums  over  the  eigenfunctions  of  the  problem  of  solving  Schrodinger’s  equa- 
tion within  the  sphere,  subject  to  the  boundary  condition  of  having  a prescribed  log- 
arithmic derivative  at  the  surface  of  the  spheres.  The  other  of  the  present  writers 
(MMS),  examining  this  formulation,  noticed  that  these  infinite  sums  represented  the 
expansions  of  certain  Green's  functions,  in  series  of  the  eigenfunctions  so  set  up,  and 
it  occurred  to  him  that  one  could  write  these  Green's  functions  directly,  thus  in  effect 
finding  explicit  expressions  for  these  infinite  sums.  'I'his  has  now  been  done,  with  a 
resulting  great  simplification  of  the  formalism  of  the  method,  but  at  the  same  time 
with  an  alternative  physical  description,  which  shows  us  exactly  the  navire  of  the  final 
wave  functions  and  energy  levels.  We  shall  first  state  this  alternative  description, 
then  give  the  resulting  mathematical  formulation. 

This  alternative  physical  description  proves  to  be  the  following.  We  start  as 
before  with  the  plane  wave  outside  the  spheres.  We  now  set  up  a solution  of  Schr6- 
dinger's  equation  within  the  sphere  for  an  arbitrary  energy  E,  soon  to  be  determined. 
For  any  such  energy  E,  we  can  at  once  set  up  a solution  which  is  continuous  with  the 
plane  wave,  but  the  derivative  will  not  in  general  be  continuous.  We  can  then  compute 
the  expectation  value  of  the  energy,  for  the  wave  composed  of  this  solution  of  Schro- 
dinger's  equation  inside  the  sphere,  and  the  plane  wave  outside.  This  expectation 
value  will  have  as  two  of  its  terms  E^  times  the  integral  of  the  square  of  the  plane 
wave  over  the  region  outside  the  spheres,  and  E times  the  integrad  of  the  square  of 
the  wave  function  over  the  inside  of  the  spheres.  There  is  a third  contribution,  how- 
ever, coming  from  the  discontinuity  of  slope  of  the  wave  function  on  the  surface  of  the 
sphere.  This  discontinuity  in  slope  can  he  considered  a limiting  form  of  an  infinite 
second  derivative,  which  integrates  to  a finite  contribution  to  the  kinetic  energy,  when 
we  integrate  over  a thin  spherical  shell  at  the  surface  of  the  atom.  We  now  take  this 
resulting  expectation  value  of  the  energy,  and  make  the  requirement  that  it  be  station- 
ary when  we  vary  the  energy  E inside  the  spheres.  When  we  do  this,  we  find  that  the 
resulting  equation  can  b«»  satisfied  in  a simple  fashion;  the  expectaxion  value  of  the 
energy  is  set  equal  E,  the  energy  for  which  we  are  solving  Schrddinger 's  equation  in- 
side the  sphere.  In  other  words,  we  have  achieved  an  exact  solution  of  SchrSdinger's 
equation  for  the  energy  E inside  the  spheres,  and  the  discrepancy  in  energy  resulting 
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from  the  fact  that  the  energy  of  the  plane  wave  part  of  the  wave  function  outside  the 
spheres  is  E^  is  compensated  by  an  equal  and  opposite  contribution  from  the  surface 
integral. 

The  mathematical  formulation  of  this  statement  of  the  method  is  very  simple, 
and  we  shall  give  its  essentials  in  a moment.  The  interesting  fact,  however,  is  that 
we  can  take  the  solutions  of  Schrodinger 's  equation  for  the  energy  E,  as  encountered 
in  this  method,  and  expand  them  in  terms  of  the  eigenfunctions  used  in  the  earlier 
method,  where  we  set  up  eigenfunctions  of  the  problem  with  a prescribed  logarithmic 
derivative  at  the  surface  of  the  sphere.  When  we  do  this,  we  can  easily  show  that  this 
method  is  identical  with  the  other  one,  the  infinite  sums  encountered  in  the  earlier 
form.ulation  becoming  formally  identical  with  much  simpler  functions  found  in  the  new 
method.  We  see,  in  other  words,  that  the  variation  problen.  set  up  in  the  older  ver- 
si<^n  of  the  method,  coupled  with  the  expansion  of  the  solution  of  Schrodinger 's  equation 
inside  the  sphere  in  terms  of  an  infinite  sum  of  eigenfunctions  all  having  a prescribed 
value  of  logarithmic  derivative  at  the  surface,  so  as  to  produce  continuity  of  slope  as 
well  as  function  with  the  plane  wave,  has  automatically  built  up  a solution  inside  the 
sphere  which  docs  not  have  continuity  of  slope  with  the  plane  wave.  It  does,  however, 
have  the  property  we  have  just  mentioned,  that  it  satis'^es  Schrodinger's  equation  ex- 
actly inside  the  spheres,  and  uses  the  discontinuity  of  slope  to  compensate  for  the  fact 
that  it  does  not  satisfy  it  exactly  outside  the  spheres.  It  still  appears  that  these  func- 
tions should  form  a very  good  starting  point  for  a secular  problem  involving  a number 
of  such  augmented  plane  waves. 

The  mathematical  formulation  is  now  very  simple-  Let  an  atomic  sphere  be 
located  at  the  origin,  with  a radius  K.  Let  U|(E;  r)  be  the  radial  part  of  a solution  of 
Schrodinger's  equation  for  the  interior  of  the  sphere,  with  aximuthal  quantum  numbe’' 
i.  and  energy  E.  and  regular  at  the  origin.  Then,  using  the  familiar  expansion  of  the 
plane  wave,  the  solution  of  Schrodinger's  equation  for  energy  E inside  the  sphere, 
continuous  with  the  plane  wave  of  unit  magnitude  outside,  can  be  written 


i u^(E;r) 

X(f)(2i  + 1)  i j/kK)  P,(cos  0)  ; (2.  i) 

u^(E;r) 


when  r = R,  this  at  once  reduces  to  the  expansion  of  the  plane  wave  at  radius  R.  We 
can  then  show  easily  that  the  contribution  of  the  surface  discontinuity  to  the  expecta- 
tion value  of  the  energy  is 


4ttR^  I(i)(2i  + 1)  j/(kR) 


d In  U|(E;r)  d In  j|(kr) 


dr 


dr 


r = R 


(2.2) 


ft 
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The  last  term  in  the  square  bracket  can  be  eliminated,  if  we  choose,  by  using  the  re- 
lation 

X(f)(2f  + 1)  J|(2)  = 0 . 


which  can  be  proved  from  the  properties  of  spherical  Bessel  functions.  The  contribu- 
tion of  the  plane  wave  to  the  discrepancy  between  the  expectation  value  of  the  energy 
and  E^  is  - E).  v/here  fi  is  the  vobimp  of  space  outside  the  sphere.  Since  our 

condition  on  E is  that  this  discrepancy,  plus  the  contribution  of  the  surface  discon- 
tinuity to  the  energy,  must  add  to  zero,  our  condition  for  the  energy  is  found  by  setting 
the  expression  (2)  equal  to  fJ(E  - E^).  We  are  left,  then,  with 


, , d In  u,(E;r) 

n(E  - Ey)  = 4ttR^  $;(f)(2f  + 1)  j/(kH)  


r = R 


(2.  3) 


as  the  condition  for  the  energy.  We  can  solve  this  by  plotting  the  right  side  as  a func- 
tion of  E,  and  finding  the  intersection  of  this  curve  with  the  straight  line  represented 
by  the  loft  side.  The  details  of  the  resulting  conclusions  regarding  the  energy  levels 
are  given  in  the  next  contribution  to  this  Progress  Report. 
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3.  ENERGY  RELATIONS  IN  THE  AUGMENTED  PLANE  WAVE  METHOD 


In  order  to  discuss  the  solutions  of  the  energy  (Eq.  2.  3 of  this  Report),  we 
plot  both  the  left-hand  and  right-hand  members  of  the  energy  expression  as  functions 
of  the  energy,  the  intersections  of  these  curves  giving  us  solutions.  The  left-hand 
member  is  simply  a straight  line  with  slope  Q which  intersecis  the  energy  axis  at 
E = E^.  The  right-hand  member  is  a superposition  of  the  curves  d In  u^(E;  r)/drj^_j^ 
taken  with  the  weight  4TrR^{2f  + 1)  j^^(kR).  Any  such  curve  will  have  infinites  for 
those  values  of  energy  which  cause  u^(E;  R)  to  equal  zero.  Since  the  slope  of  such  a 
curve  must  be  everywhere  negative,  as  we  can  show,  the  curve  will  have  a strong  re- 
semblance to  a cotangent  curve.  This  behavior  is  consistent  with  what  we  should  ex- 
pect if  wc  approximated  our  wave  functions  U|  by  the  WKB  method.  For  energies 
which  are  very  negative  the  solutions  to  the  Schrouingei  equation  are  essentially  a 
combination  of  the  exponentials 

^-yiiir  +yiiir 

£ and  . 

r r 

Thus  for  r = R the  wave  function  is  practically  the  increasing  exponential,  and  our 
curve  is  essentially  the  parabola  - l/R.  As  we  come  to  the  first  bound  state 

the  increasing  exponential  comes  in  with  a smaller  and  smaller  contribution  until  it 
disappears  completely,  and  then,  soon  after,  the  wave  function  vanishes.  For  very 
negative  energies  the  bound  state  occurs  at  the  intersection  of  the  curve  and  the  parab- 
ola - -v/lE  I - l/R.  This  behavior  repeats  for  the  other  very  tightly  bound  states.  As 
the  energy  increases  further  the  decreasing  exponential  becomes  increasingly  impor- 
tant at  r = R,  so  that  between  two  infinities  in  the  curve,  the  slope  changes  more 
gradually.  As  the  energy  increases  still  further,  the  many  oscillations  in  the  wave 
functions  reach  the  point  r = R,  and  the  curve  takes  on  more  and  more  a cotangent 
behavior.  In  the  region  about  zero  energy  the  shape  of  the  curves  will  depend  strongly 
on  the  shape  of  the  potential.  In  the  region  of  large  positive  energy  the  wave  function 
will  essentially  be  a spherical  Bessel  function  with  a small  admixture  of  spherical 
Neumann  function,  and  the  curve  is  sensibly  a cotangent  curve.  This  description  is 
summarized  in  Fig.  3-1.  As  the  angular  momentum  increases,  the  wave  function 
must  go  out  further  and  further  in  distance  before  it  builds  up.  It  also  dies  out  more 
gradually  and  so,  for  a bound  state,  can  tie  on  only  to  a decreasing  exponential  which 
dies  out  more  gradually.  Collecting  the  consequences  of  raising  the  angular  momen- 
tum we  can  say:  the  first  bound  state  occurs  at  higher  energy  the  higher  the  angular 
momentum:  the  effect  of  the  decreasing  exponential  is  felt  more  and  more  even  in 
lower  energies;  the  oscillations  in  the  wave  function  are  felt  at  r = R earlier  in  en- 
ergy; and  the  shape  of  the  curve  about  zero  energy  depends  less  and  less  upon  the 
shape  of  the  potential  in  the  inner  regions.  For  very  high  angular  momenta  there  are 
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u!/E;R.) 

— I'B  D \ ^ function  of  E for  small  1 

u./E;R.) 


no  bound  states  and  for  negative  energies  the  curve  is  essentially  equal  to  the  angular 
momentum  quantum  number,  f,  until  it  drops  gradually  to  zero  and  then  down  to  minus 
infinity,  becoming  for  all  higher  energies  a cotangent  curve.  This  occurs  because  for 
high  energies  the  wave  function  is  essentially  a spherical  Bessel  funcuon  whose  be- 
havior is  just  as  stated. 

We  can  now  picture  the  curve  of  energy  of  the  wave  function  as  a function  of 

the  energy  of  the  plane  wave.  This  curve  is  traced  out  if  we  plot  the  energy  values  at 

which  the  straight  line  S1(E  - E cuts  the  cotangent -like  curve  4jtR^  ^(i)(2f  + 1)  j^^(kR) 

d In  u^(E;r)/drj^  - R energy  of  the  plane  wave  increases  indefinitely.  If  we  rail 

the  continuous  portion  of  the  curve  between  two  singularities  (two  adjacent  singularities 

will  in  general  belong  to  wave  functions  of  different  angular  momenta)  a "branch  of 

the  curve,  then  to  each  such  branch  of  the  cotangent -like  curve  there  corresponds  a 

branch  of  the  E vs.  E curve.  Each  of  these  branches  falls  between  lines  of  constant 

o 

energy  E,  and  these  lines  occur  at  all  those  energies  for  which  any  one  of  the  wave 
functions  Uj(E;  R)  vanishes.  As  the  plane  wave  increases  in  energy  the  cotangent -like 
curve  changes  shape  but  each  branch  remains  confined  between  the  two  asymptotes 
which  bound  it  and  these  do  not  change  position  as  the  plane  wave  increases  in  energy. 
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At  the  same  time  the  straight  line  I2(E  - E^)  moves  down  the  energy  axis  with  fixed 
slope  cutting  the  axis  at  E = E^.  Let  us  picture  what  happens  in  a typical  branch  cen- 
tered about  E = E'.  When  E^  is  very  much  less  than  E'  the  straight  line  will  intersect 
the  branch  very  close  to  its  left-hand  asymptote  because  the  straight  line  will  have 
achieved  quite  a laige  value  before  it  intersects  this  branch,  and  only  those  portions 
of  the  branch  near  the  asymptote  can  reach  that  high.  As  E^  approaches  E',  the  straight 
line  will  cut  the  branch  at  points  which  arc  closer  to  the  right-hand  asymptote.  How- 
ever, suppose  Ey  nears  a value  which  causes  j^(kR)  to  vanish  and  I is  the  angular  mo- 
mentum quantum  number  equal  to  the  index  of  the  wave  function  u^  which  generates  the 
left  asymptote.  Then  as  the  plane  wave  energy  nears  that  point,  the  contribution  of 
d In  Uj(E;  r)  dr|^  _ to  the  cotangent -like  curve  nears  zero  except  at  the  points  where  it 
has  iis  singularities.  At  these  points  the  branches  to  cither  side  of  the  singularity  will 
be  slowly  changing  except  very  close  to  the  asymptote  whore  they  tend  to  flatten  out 
against  it  and  then  sharply  achieve  singular  value.  The  closer  the  plane  wave  energy  to 
the  point  which  causes  this  behavior  the  more  pronounced  the  behavior  becomes.  If 
this  point  at  which  u^(E;  R)  vanishes  is  sufficiently  distant  from  the  brar.ch,  the  intersec- 
tion of  the  branch  with  the  straight  line  nears  the  energy  of  the  asymptote  as  the  plane 
wave  energy  increases.  This  causes  the  E vs.  E^  curve  to  dip  down  to  its  asymptote 
and,  as  goes  beyond  this  point  the  curve,  rises  again  to  resume  its  regular  behavior. 

The  result  is  that  the  E vs.  curve  is  tangent  to  its  lower  asymptote  at  such  a point. 
Now  let  us  suppose  that  the  plane  wave  energy  has  passed  such  a point  and  is  nearing 
E',  the  location  of  the  branch.  Then  the  intersection  moves  across  the  branch  until, 
when  E^  is  very  much  greater  than  E',  the  intersection  is  made  with  values  of  energy 
near  the  right-hand  asymptote.  The  behavior  of  the  energy  curve  when  E^  yields  a 
root  of  a spherical  Bessel  function  with  E^  » E'  the  same  as  that  in  the  previous  case 
mentioned  except  that  now  the  energy  curve  dips  up  to  meet  the  upper  asymptote  of  the 
energy  branch. 

Thus  if  we  draw  lines  of  constant  plane  wave  energy  at  all  of  those  energies 
for  which  any  one  of  spherical  Bessel  functions  vanish,  then  when  these  lines  intersect 
the  asymptotic  lines,  the  energy  curve  will  be  tangent  to  the  asymptotic  line  at  this  in- 
tersection only  if  both  these  lines  correspond  to  the  same  angular  momentum  and  also 
if  this  intersection  is  sufficieiatly  distant  from  the  region  in  which  E = E^ . Every 
spherical  Bessel  function  except  that  of  order  one  vanishes  at  E^  = 0.  Thus,  every 
energy  curve  will  start  off  at  E^  = 0 touching  its  lower  asymptote  if  that  asymptote 
corresponds  to  any  angular  momentum  other  than  zero  and  also  if  that  asymptote  is 
Sufficiently  distant  from  zero  energy.  There  is  a ramification  of  this  peculiar  event  in 
the  properties  of  the  wave  function  which  shall  not  be  examined  here.  We  now  discuss 
the  shape  of  the  curve  for  negative  energies. 

It  is  clear  that  the  branches  with  very  negative  energies  will  intersect  the 
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Fig.  3-2 

Energy  of  an  augmented  plane  wave  as  a function 
of  the  energy  of  its  plane  wave  portions.  (Circled 
points  denote  points  of  tangency.  ) 

straight  line  ve 'y  near  the  right-hand  asymptote  for  all  plane  wave  energies  (we  dis- 
count for  the  moment  the  peculiar  behavior  mentioned  above).  For  the  less  tightly 
bound  states,  however,  there  remains  the  possibility  that  if  the  branch  is  sufficiently 
depr»*ssed  the  energy  curve  can  start  off  near  the  left  asymptote  and  then  rise  to  the 
right-hand  asymptote  very  n;uch  as  the  positive  energy  branches  dc.  At  least,  there 
can  be  a remnant  of  such  behavior  for  these  branches.  We  have  already  mentioned 
that  the  asymptotes  for  very  negative  energies  are  very  close  to  the  energies  of  the 
strongly  bound  states  so  we  see  that  the  bound  state  energies  are  essentially  independ- 
ent of  E^.  It  is  also  clear  that  the  peculiar  behavior  associated  with  the  vanishing  of  a 
spherical  Bessel  function  will  occur  also  for  negative  energies,  except  perhaps  for 
those  states  near  zero  energy.  The  energy  curve  is  shown  in  Fig.  3-2.  The  shape  of 
this  curve  along  the  line  E = E^  is  reminiscent  enough  of  energy  bands  and  the  true 
periodic  problem  to  suggest  that  these  portions  of  the  energy  curves  are  rough  approxi- 
mations to  a band  structure.  It  also  suggests  which  of  the  augmented  plane  waves  to 
choose  when  we  wish  initial  functions  for  our  periodic  problem  when  we  attempt  solu- 
tions for  the  higher  bands. 

Since  the  straight  line  hss  the  slope  Q,  it  should  be  clear  that  the  effect  of 
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increasing  SI  is  to  make  the  onset  of  the  rise  of  the  energy  curve  near  E = E^  more 
sudden  and  to  make  he  curve  everywhere  else  practically  identicad  with  the  asymptotic 
values.  The  effect  for  negative  energies  is  a tendency  to  remove  the  structure  from 
the  energy  curves  around  the  region  k = 0.  However,  it  should  be  noted  that  the  effects 
o'”  making  P.  larger  is  slightly  compensated  in  the  case  of  negative  energy  by  the  de- 
crease in  R,  the  compensation  arising  from  a depression  of  the  cotangent -like  curve 
in  that  energy  region.  There  is  still  a third  factor.  The  curve  will  take  on  positive 
energy  behavior  much  later  in  negative  energy.  This  effect  also  tends  to  compensate 
the  increase  in  P. 


M.  M.  Saffren 
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4.  ENERGY  BANDS  IN  CHROMIUM 


The  calculation  of  matrix  elements  for  the  augmented  orthogonalized  plane 
wave  method  is  in  progress.  In  order  to  calculate  kinetic  energy  integrals  and  Fourier 
transforms  of  the  atomic  orbitals  it  is  most  convenient  to  express  these  atomic  or- 
bitals as  simple  analytic  functions.  In  conjunction  with  A.  J.  Freeman,  the  lon-core 
atomic  orbitals  of  chromium  have  been  so  expressed.  At  present  this  process  is  being 
carried  out  for  the  modified  3d  orbital  discussed  in  the  last  Progress  Report. 

R.  H.  Parmenter 
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5.  A REVISED  VALENCE  BOND  THEOHY 


Recent  work^^^  has  shown  that  it  is  possible  t^  carry  through  completely  non- 
empirical  caloulations  on  simple  molecules  using  a revised  form  of  the  valence  bond 
(VB)  theory;  since  this  method  is  to  be  developed  during  the  coming  year  it  seems  de- 
sirable to  repor'  br'efly  o,i  what  has  already  been  accomplislied  and  on  plans  for  fur- 
th.er  work. 

The  conventional  VB  theory  proceeds  in  the  following  way:  electrons  are  alotted 
to  the  AO's  a,  b,  c.  --  in  some  way  (usually  there  is  just  one  electron  in  each  orbital 
but  in  so-called  "polar"  functions  doubly  filled  (-  ve)  and  empty  (+  ve)  orbitals  occur 
--  these  are  not,  however,  invoked  in  the  usual  approximation).  Spin  factors  are  then 
added  to  give  a product  function  such  as  U = a(l)  a(l)  b(2)  p(2)  c(3)  p(3)  . . . and  an 
antisymmetrized  product  (AP)  is  constructed;  <{)  = j/u  = |(aa)(bp)  ...  I say.  A large 
number  of  eligible  4>'s  arises  because  there  are  many  possible  spin  allocations;  but 
these  4>'s  may  be  grouped  into  combinations  which  are  spin  eigenfunctions  (SE's)  and, 
in  constructing  an  approximate  wave  function,  only  those  SE's  with  similar  eigenvalues 
of  the  operators  S and  S^  can  be  mixed,  it  is  usual  to  consider  a singlet  ground  state 
and  to  employ  SE's  formed  by  coupling  the  spins  in  pairs  --  these  arc  the  "structures" 
(we  shall  denote  them  by  4i's)  of  VB  theory.  Finally,  approximate  wave  functions 
(J's)  are  constructed  in  the  basis  of  4''s.  By  analogy  with  the  Heitler-London  calcula- 
tion on  H^.  a single  structure  might  be  expected  to  describe  a system  of  "bonds”,  each 

bond  involving  orbitals  with  paired  spins:  but,  as  Slater  has  pointed  out  for  the  IF 

(3)  ^ 

case,  ' spin-pairing  alone  is  no  criterion  for  bonding,  for  a change  in  character  of 

the  AO's  may  lead  to  a v|i  predicting  strong  repulsion.  A VB  energy  calculation  gen- 
erally becomes  feasible  only  with  the  neglect  of  overlap  between  different  AO's;  in 
other  words,  conventional  energy  calculations,  employing  a basis  of  VB  structures, 
can  only  be  rigorously  validated  by  a preliminary  orthogonalization  of  all  orbitals;  and 
our  experience  with  the  hydrogen  molecule  shows  Just  how  crucial  a step  this  may  be. 

In  the  revised  theory  we  do  indeed  start  from  a set  of  localized  but  orthogonal 

orbitals  a,  5,  c,  --  (we  shall  call  these  AO's),  constructed  from  the  basic  AO's  a,  b, 

(4) 

c,  --  by  the  method  of  Lowdin.  It  is  then  possible  to  preserve  the  essential  formal- 
ism of  the  VB  approach,  while  completely  avoiding  the  approximations  of  the  usual 
treatment.  But  the  characteristic  features  of  the  VB  method  are  found  to  have  com- 
pletely disappeai  ed:  the  "VC"  structures  of  the  new  theory  (i.  c.  , structures  defined 
over  AO's  rather  than  AO’s)  do  not  describe  systems  of  chemical  bonds  in  any  sense 
whatever;  in  fact,  the  conclusions  reached  by  Slater  in  the  case  of  are  found  to  be 
quite  general  --no  single  "covalenf  VB  structure  can  describe  chemical  binding,  pre- 
dicting instead  strong  repulsion  at  adl  internuclear  distances.  Moreover,  the  situation 
cannot  be  restored  by  the  introduction  of  any  number  of  alternative  non-polar  VB  struc- 
tures: in  order  to  describe  a bond  involving  orbitals  a and  D it  is  absolutely  essential 
to  admit  polar  VB  structures  so  that  the  final  wave  function  will  contain,  with  substan- 
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tial  weight,  pairs  of  structures  which  differ  by  a "shift  of  charge"  between  the  orbitals 
a and  H.  The  bond  a— b therefore  arises  as  an  inter-configurational  effect  from  the 
pairs  of  structures  associated  pictorLally  with  a charge  "hop"  between  atoms  a and  b. 
This  conclusion  is  readily  interpreted  in  terms  of  the  molecular  charge  density  func- 
tion. The  basic  AO's  have  been  replaced  by  distorted  AO's  (AO's)  and,  m the  absence 

—2  ~2T  — —Z 

of  polar  structures,  J can  only  lead  to  a density  of  the  form  P = a + E + c + ...; 

but  the  AO's  are  clearly  so  distorted  as  to  give  a particularly  low  charge  density  in 
the  bond  regions.  'I'he  high  internuclear  charge  density,  with  which  binding  is  asso- 
ciated, is  introduced  by  the  cross  terms  in  $ ^ which  come  from  structures  related 
by  a charge  hop. 

Calculations  by  the  "VB"  method,  obviously  require  the  use  of  both  non-polar 
and  polar  (possibly  including  multiply -polar)  VB  structures  and  consequently  the  cal- 
culation of  matrix  elements  for  all  types  of  pair.  This  problem  has  been  solved  by  an 

(5) 

extension  of  the  cycle  method  of  Eyring  and  Kimball;  the  results  will  be  presented 
here  without  proof. 

Any  singlet  structure  is  a combination  of  ^'s  which  shows  antisymmetry  in  the 
spin  factors  of  various  orbital -pairs;  in  the  Rumer  notation  such  a structure  might  be 

a b 
h c 

t * 

g d 

f e 

where  the  pairs  arc  shown  linked  and  the  arrow  indicates  the  order  of  the  spin  factors 
(a  -►  p)  in  the  leading  (or  "parent")  4>  of  the  SE  combination.  * Triplet,  quintuplet,  -- 
structures  arise  when  antisymmetry  exists  for  one,  two,  --  less  than  the  full  number 
of  possible  pairs,  the  combination  being  symmetric  in  the  spins  of  all  the  unpaired 
orbitals.  But  (at  least  for  not  too  high  multiplicity)  such  an  SE  is  expressible  in  terms 
of  a small  number  of  completely  paired  structures  in  which  some  pairs  are  now  as- 
sociated with  spin  symmetry.  Thus,  the  quintuplet  SE 
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Reversal  of  an  arrow  simply  changes  the  sign  of  a iV:  directions  are  therefore  unim- 
portant except  in  determining  the  sign  of  a matrix  element. 
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is  equivalent  to  the  combination 


where  we  now  distinguish  the  two  symmetry  types  by  + and  New  configurations 
may  be  introduced  by  substituting  a new  orbital  for  one  of  those  in  the  standard  set 
a,  b,  c,  --  (indicated,  for  example,  by  f(k))t  in  particular,  polar  structures  occur  when 
the  orbitals  at  the  ends  of  a (-)  link  are  identified.  Thus  a b(a)  indicates  that  or- 

bital b has  been  replaced  throughout  by  a - giving,  when  the  orbitals  are  differently 

localized,  a b^;  but  the  possible  spin  interchange  between  linked  orbitals  becomes 

redundant  in  this  case  and  the  link  should  be  distinguished  as  a "dummy"  or  (o)  link. 

The  matrix  element  for  any  pair  of  spin-paired  structures  and  the  full  Hami’tonian  H 
may  now  be  written  down  by  inspection  of  a cycle  diagram  (set  up  by  the  method  of 
Eyring  and  Kimball,  disregarding  new  link -types  and  orbital  substitutions),  all  links 
being  inserted  so  as  to  give  a continuous  non-crossing  chain  in  each  cycle.  Thus,  for 
two  typical  triplet  state  structures  with  the  following  pairing  schemes 

iVj:  a-^b(a)  c-^d(c)  e — f g-^h(g)  i-— j k-^1  m— n (iriply  polar) 

4<2t  a— b(a)  c-^f  d-^e  h-— j g~i  l-^m  k— ^n  (singly  polar) 

the  cycle  diagram  is 


where,  for  instance,  d(cd)  means  that  d is  substituted  by  c in  but  left  unchanged 
in  4<2-  noticed  that  vertical  links  refer  to  and  diagonal  to  ii<2-  ^he  re- 

quired matrix  element  is  then 

(4»,|HI4i,)  = [q*  - Vx.  K'  - Vk  K’  1 

1 L ^ pq  pq  ^*^pq  pqj 

where  O = (abed  --  IHI  abed  --),  K = (pqfgiqp)  (g  = l/r  ,),  n is  the  number  of  cycles, 

pq  1 ‘ 

nj  being  the  number  containing  (O)  links  only,  and 
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p,  q in  different  cycles  p,  q in  the  same  cycle 


pq 


= 1/2 


pq 


= 0 


pq 


pq 


pq 


pq 


= 1 p,  q same  side. 

= 0 p,  q different  sides. 

= t 1 according  as  the  number 
of  (-)  links  in  each  path 
between  p and  q is  even  or 
odd  --  (0)  links  regarded 
as  equivalent  to  (-). 

= 0 p,  q joined  by  a (0)  link  or 
by  a chain  of  (O)  links  only. 


There  is  no  no  cunii  ibution  if  any  cycle  contains  an  odd  number  of  (+)  links.  I'he 
primes  on  the  Q's  and  the  K's  indicate  that  the  standard  set  of  orbitals  a,  b,  c,  --  may 
be  modified,  in  either  or  or  in  both,  by  orbital  substitutions!  the  sums  extend 

over  distinct  pairs,  p,  q and  vanishes  unless  the  pair  p,  q includes  the  one  or  two 

orbitals  in  which  and  4*2  oaay  differ.  Application  of  the  rule  is  very  simple;  in  the 
present  case  we  see  immediately  (4,j|Hl4<2)  = 2^  [Qbdh(aa,  cd,  gh)  - ^ K^glcd,  gh)] 
where  Qb^jh(aa,  cd,  gh)  = (abcdcfgR  --iHlaSc^efgli  --)  and  K^glcd,  gh)  = (Sflgl^^), 
the  raised  letters  indicating  the  substitutions  to  be  made.  These  quantities  are  always 
immediately  expressible  in  terms  of  one  - anu  two-particle  integrals;  in  the  present 
case  Q'  reduces  simply  to  (cgigidh)  --  a generalized  Coulombic  integral  which  might, 
in  the  present  notation,  have  been  called  *^(jb^cd,  gh). 

The  method  outlined  is  quite  general;  but  in  dealing  with  the  VB  structures, 
where  the  orbitals  which  we  have  labelled  a,  b,  c,  --  are  in  fact  the  AO's  a,  F,  c,  -- 
certain  features  of  the  subsequent  energy  calculations  are  particularly  noteworthy. 
They  are 

(1)  The  only  large  integrals,  which  determine  the  essential  structure  of  the 

___ 

final  secular  equations,  are  of  the  types  = (alhla),  = (alhib)  and  = 

(abigiab),  a and  P being  the  two  integrals  which  would  dominate  a molecular  orbital 
treatment.  Ail  generalized  Coulombic  integrals  (such  as  J^^(cd,  gh),  above)  and  all 
exchange  and  generalized  exchange  integrals  are  very  much  smaller,  when  defined 
over  AO's,  than  in  conventional  treatments. 

(2)  Intra-configurational  matrix  elements  contain  only  a's  and  J's  (plus  insig- 
nificant exchange  terms);  they  arise  from  groups  of  structures  which  describe  non- 
bonded  atoms  interacting  as  independent  units  --  Q'  = (abbd  - IHlabbd  --),  for  in- 
stance, represents  the  energy  of  a system  with  one  electron  on  a,  two  on  b,  none  on 
c and  one  on  d,  including  mutual  repulsions. 


h is  the  "effective"  Hamiltonian  for  one  electron  in  the  field  of  nuclei  and  closed  shells. 
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(3)  Inter -configurational  elements  contain  only  ]5's  (plus  insignificant  exchange 
terms),  arising  from  pairs  of  structures  differing  by  a shiit  of  charge  between  the 
atoms  about  which  a and  H are  localized,  and  being  simply  the  "resonance"  or  "hopping" 
integral  of  MO  aieory.  It  is  the  presence  of  these  elements  alone  which  allows  us  to 
account  for  chemical  binding,  within  the  VB scheme. 

Detailed  numerical  calculations  are  now  in  progress  on  the  it  electron  systems 
of  three  organic  molecules  --  ethylene,  cyclo-butadience  and  benzene.  The  first 
case  is  trivial;  but  is  useful  because  it  admits  a comparison  between  non-empirical 
estimates  of  a and  jJ  and  the  best  empirical  values  (from  ionization  potential  and 
spectral  data).  The  second  case  is  important  because  it  makes  possible  a compari- 
son of  the  VB procedure  with  the  complete  MO  treatment  given  by  Craig,  who  used 
for  the  singlet  states,  for  instance,  a basis  of  twenty  MO  configurations  (in  our  ap- 
proach we  employ  twenty  VB  structures).  In  the  third  case  it  will  be  possible  to  make 

a similar  comp.arison  with  the  (much  heavier)  benzene  molecule  calculations  of  Parr, 

(7) 

Craig  and  Ross,  ' ' who  employ  limited  configuration  interaction  (using,  for  instance, 
the  45  lowest  singlet  configurations);  it  is  felt  that  their  results  may  be  more  easily 
reached  and  extended  by  using  the  VB  approach.  In  all  these  cases  it  is  possible  to 
rigorously  compare,  for  the  first  time,  the  rapidity  of  convergence  of  two  approxima- 
tion techniques  which  proceed  from  completely  different  starting  points. 

Results  of  some  importance  have  already  emerged.  In  both  ethylene  and  cyclo- 
butadiene it  is  instructive  to  set  up  suitable  approximate  wave  functions  in  terms  of 
the  conventional  VB  structures  (such  functions  being  here  determined  completely  by 
symmetry):  these  cannot  be  properly  handled  in  the  usual  treatment  --  but  the  Heitler- 
London  calculation  suggests  that  they  are  plausible  wave  functions;  and  we  may  now 
expand  them  in  terms  of  VB  structures  and  carry  out  a rigorous  energy  calculation. 

It  turns  out  that  the  functions  constructed  in  this  way  resemble  quite  closely  the  best 
approximations  which  can  be  obtained  in  the  given  basis;  they  are  considerably  su- 
perior to  corresponding  one -configuration  MO  functions.  The  expanded  VB  functions 
somewhat  underestimate  the  weights  of  the  polar  VB  structures,  almost  completely 
omitting  the  small  multiply-polar  contributions  which  feature  in  the  best  approxima- 
tions. On  the  other  hand,  the  (one -configuration)  MO  method  grossly  ovei  weights  all 
the  polar  structures  and  strong  configuration  interaction  is  necessary  to  correct  this 
obvious  defect.  It  seems  likely  that  these  results  may  be  general  and  that  the  con- 
ventional VB  structures  may  be  of  great  value  in  indicating  those  groupings  of  VB 
structures  which  are  most  useful  as  "units"  from  which  good  approximate  wave  func- 
tions may  be  rapidly  constructed;  a preliminary  energy  calculation  could  then  be 


♦The  "separation"  of  w electrons  from  those  of  the  closed  shells  is  a valid  first  ap- 
proximation provided  the  effective  Hamiltonian  h is  properly  defined. 
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made  in  the  "VB  basis"  (using,  of  course,  the  VB  methods  we  have  developed)  and  sub- 
sequent refinement,  where  practicable,  could  be  made  in  the  full  basis  of  VB  structures. 
This  final  step  would  usually  be  precluded  by  the  large  number  of  structures  to  be  con- 
sidered; for  this  reason  the  calculations  on  benzene,  which  should  give  more  conclu- 
sive evidence  of  the  efficacy  of  the  simpler  procedure,  are  of  special  interest. 
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6.  ANTIFERROMAGNETISM 


In  the  preceding  Progress  Report  an  investigation  of  the  Mn^—  — Mn^^  sys- 

tem was  presented.  As  stated  before,  the  object  of  this  work  was  to  find  out  how  the 
spins  of  two  magnetic  ions,  whose  wave  functions  do  not  overlap,  arc  coupled  when 
they  are  separated  by  an  intervening  non-magnetic  ion.  Within  the  framework  of  the 
approximation  made  there,  it  was  found  that  the  Mn  ions  coupled  fcrromagnctically, 
i.  e.  , the  net  spins  of  these  ions  were  found  to  be  parallel  in  the  lowest  state  of  the 

In  order  to  examine  this  problem  of  spin  coupling  in  grea'.er  detail  a much  sim- 
pler model  was  investigated.  This  simplified  system  is  shown  in  Fig.  6-1.  The  one- 


Uj  s,p  u^ 

He^  He  He^ 

•^R^  =2.5  a.  u. 


Fig.  6-1 

electron  orbitals  used  were  gaussians  localized  about  the  various  centers.  Let  u.  de- 

4.  * 

note  an  s-like  gaussian  about  the  left-hand  He  and  u.  an  s-like  gaussian  about  the 
right-hand  He  . Let  s and  p denote  an  s-like  and  p-like  gaussian  respectively  about 
the  central  He.  Two  types  of  interactions  were  considered  and  their  magnitudes  com- 
pared. The  first  type  was  the  spin  coupling  of  the  electrons  on  the  two  end  He^  ions 
brought  about  by  configurations  in  which  the  electrons  on  the  central  He  were  excited 
into  higher  states  of  that  atom  but  were  not  transferred  to  one  of  the  end  He^  ions. 
These  configurations  are  listed  in  Table  6-1. 

T able  6-1 

Uj  s p u^ 

1 2 0 1 .1,^^ 

1 1 1 1 +sp 

1 0 2 1 +PP 

The  second  type  of  interaction  investigated  was  the  spin  coupling  of  the  two  end  He^ 
ions  brought  about  by  the  superexchange  configurations,  in  which  one  electroii  is 
taken  from  the  central  He  and  placed  on  one  of  the  end  He^  ions.  These  configurations 
are  listed  in  Table  6-2  where  only  s orbitals  were  considered  in  this  case. 

It  was  felt  th.at  by  using  non-orthogonal  orbitals  for  this  simple  model  some 
light  might  be  thrown  on  the  question  of  whether  or  not  the  use  ol  orthogonal  one- 
electron  orbitals  and  a rather  limited  configuration  interaction  for  the  Mn-*-0^ — Mn^^ 
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Z 

! 


Table  6-2 

+0 


problem  was  responsible  for  the  result  found  there.  Also  the  results  for  this  simple 
system  can  be  re-expressed  in  terms  of  orthogonal  orbitals  showing  what  configura- 
tions of  orthogonal  orbitals  would  give  an  equivalent  result. 

The  overlap  matrix,  m..  = / u.*u.  dv  is  given  in  Table  6-3.  The  orbitals  u, 
and  u^  were  assumed  not  to  overlap,  i.  e.  , Uj(x)  u^(x)  = 0 for  all  x.  The  non- 
orthogonality was  carried  through  completely  no  overlap  integrals  or  three -center 
integrals  being  neglected. 


Table  6-3 


Uj 

P 

^^4 

“l 

1 6 

UjS 

0 

s 

6 1 
UfS 

0 

^4S 

p 

6 r,  0 

UjP 

1 

6 

U4P 

“4 

0 6 ^ 
U4S 

^u^p 

1 

The  analytic  form  of  the  orbitals  used  was 

u,  = 

i IT 

" Ra)^ 

- Ra)^ 

p = ze~ 

2 

r 

All  of  the  gaussian  half-width  factors  were  chosen  as  unity  in  order  to  simplify  the. 
numerical  work.  The  values  of  the  overlap  integrals  were 


6 =0. 043940 

UjS 

^u^p  = 0. 109850 
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'I'he  creation  of  the  various  states  v.hich  entered  the  calculation,  their  normalization, 
and  the  matrix  components  between  them  was  treated  by  the  spin  operator  method  dis- 
cussed in  the  last  two  Progress  Reports.  Although  this  technique  makes  the  problem 
of  non-orthogonality  easier  to  cope  with,  there  is  a considerable  degree  of  complexity 
in  the  calculation. 

Singlet  and  triplet  states  corresponding  to  each  type  of  configuration  were  set 
up  for  both  types  of  interaction  and  the  inversion  symmetry  of  the  system  was  taken 
account  of.  The  results  of  the  configuration  interaction  calculation  for  those  states 
listed  in  Table  6-1  are  given  below.  * 


Singlets 


Triplets 


H 

ss,  ss 

= 

- 

15. 574729 

II 

ss,  ss 

- 

15. 574655 

= 

- 

12. 316543 

II 

sp,  sp 

= 

- 

12. 273200 

H 

PP.  PP 

= 

- 

7. 866540 

II 

PP.  PP 

= 

- 

7. 865709 

H 

ss,  sp 

= 

- 

0. 197088 

11 

ss,  sp 

= 

4 

0. 1 14860 

H 

ss,  PP 

= 

•f 

0. 379815 

II 

ss.  PP 

4 

0. 379064 

H 

sp,  PP 

= 

0.  160435 

H 

sp.  PP 

= 

- 

0. 091609 

ss’  ^sp' 

= 

+ 

0. 01  1763 

<*''ss’  '^sp) 

s 

- 

0. 006861 

1 , + ) 
ss  PP 

= 

- 

0. 011883 

(+  > + ) 
^S5  pp 

= 

0. 000000 

sp  ^pp' 

= 

- 

0. 000047 

sp  ^pp' 

= 

4- 

0. 006790 

The  terms  ('I'jj.  denote  the  overlap  integrals  between  the  various  configurations. 
The  lowest  root  of  the  secular  equation  for  the  two  multiplicities  was 


E(S  = U)  = - 15.  593377 
E(S  = 1)  rr  - 15.  593265 

The  corresponding  eigenstates  are 

^.{.S  = 0)  = 0,998‘>30-l-  + 0.003792^1/  - 0.048995  ^ 

^ ss  ^sp  ^pp 

4/{S  = 1)  = 0.998549  + - 0.002158  + - 0.048979  + 

^ ss  ^sp  ^pp 


^Atomic  units  are  used  throughout. 
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It  is  apparent  that  although  the  singlet  state  lies  lower  than  the  triplet,  the  two  eigen- 
states i'.re  very  r.early  degenerate.  One  can  conclude,  at  least  for  this  model,  that  the 
configurations  given  in  lable  6- I have  very  little  to  do  with  the  energy  dependence  of 
the  spin  orientations. 

'I'he  results  for  the  superexchange  configurations  listed  in  Table  6 -2  are  given 

below. 


Singlets 

T riplct 

s 

00 

= - 15.  574729 

**00 

= - 15. 

574655 

01 

= + 1.  181828 

»01 

= - 1. 

181015 

1 1 

= - 15.  412096 

"ll 

= - 15. 

377591 

The  overlap  integral  between  the  two  superexchange  configurations  was  the  same  for 
both  singlet  and  triplet  cases  and  was  found  to  be 

(4-0.  4',)  = - 0.  062140 

The  lowest  root  of  the  secular  equation  for  each  multiplicity  was 

E(S  = 0)  = - 15. 714038 
i:(S  = 1)  = - 15.  703501 

The  corresponding  eigenstates  were 

4;(S  = 0)  = 0.  804649  - 0.  545851 

+(S  = 1)  = 0.  824139  - 0.  517483 

Clearly  the  superexchange  configurations  lead  to  a description  of  the  system 
in  which  the  singlet  or  non-magnetic  state  lies  lower  than  the  triplet  or  magnetic  state. 
A good  treatment  of  the  model  can  be  obtained  by  leaving  out  entirely  those  configura- 
tions listed  in  Table  6-1.  ii  is  interesting  to  note  the  small  amount  of  energy  required 

to  transfer  an  electron  from  the  central  He  to  one  of  the  end  He^  ions.  As  a conse- 

quence of  this  the  final  ground  state  is  almost  an  equal  mixture  of  the  two  states  and 
of  Table  6-2. 

The  elements  of  the  secular  equation  for  the  singlet  and  triplet  states  of  the 
superexchange  configurations  are  nearly  identical  except  for  the  diagonal  energy  of  the 
4<j  state  for  the  two  multiplicities.  Apparently  the  final  singlet  state  is  lower  lying  than 
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the  final  triplet  state  chiefly  because  the  adjacent  unpaired  spins  appearing  in  the 
configuration  prefer  to  set  themselves  antiparallel  rather  than  parallel.  This  effect 
cannot  be  described  solely  n terms  of  sign  of  the  Heisenberg  exchange  integral. 

When  the  approach  of  non-orthogt>nal  orbitals  is  used,  one  is  often  tempted  to 
neglect  terms  involving  the  square  or  higher  powers  of  the  overlap  integral  and  also 
tl  r ee-  and  four -center  integrals.  Such  an  aporoximation  is  thought  to  bt>  reli-abie  iii 
those  situatioiiS  where  the  overlap  integrals  are  small  as  is  the  case  for  this  problem 
since 


J,  = 0.  043940 

UjS 

this  being  the  only  overlap  integral  appearing  in  the  superexchange  configurations. 
For  the  purpose  i f comparison  the  results  obtained  by  neglecting  these  overlap  and 
three -center  integral  terms  are  given  below. 


Singlets 


Triplets 


H 


00 


H 


01 


= + 


H. 


T 1 

(V+i)  = 


15.  542347 
1.  181454 
15.  445998 
0.  062140 


HOO  = - *5.542347 
Hqj  = + 1.  201236 
H,  j = - 15.  346550 
(4.q,4.j)  = - 0.062140 


The  lowest  root  of  the  secular  equation  for  the  singlet  states  of  the  superexchange  con- 
figurations using  this  approximation  is 


E(S  = 0)  = - 15.  705045 


This  is  to  be  compared  with  the  energy  obtained  by  carrying  the  non -orthogonality 
through  properly;  this  energy  being 

E(S  = 0)  = - 15.  714038 

The  lowest  root  of  the  secular  equation  for  the  triplet  states  obtained  by  the  above 
approximation  is 


E(S  = 1)  = - 13. 690718 

which  is  to  be  compared  with  the  rigorous  value  of 
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E(S  = 1)  = - 15.  703501 

The  approximate  singlet  triplet  energy  separation  is 

E(S  = 0)  - F,(S  = 1)  = - 0.014327 
as  compared  with  the  rigorous  separation  of 

E(S  = 0)  - E(S  = 1)  = - 0.  010537 

V\ie  see  that  even  with  such  a small  overlap  integral  the  separation  between  the  final 
singlet  and  triplet  energies,  which  is  of  primary  interest  here,  is  given  rather  differ- 
ently by  the  two  methods. 

There  is  another  interesting  piece  of  information  which  can  be  obtained  from 
this  calculation.  Suppose  we  were  tc  treat  our  non-orthogonal  orbitals  as  though  they 
were  orthogonal  and  find  the  expressions  for  the  various  matrix  elements  of  the  Hamil- 
tonain  on  that  basis.  However,  in  order  to  account  for  the  fact  that  the  orbitals  are 
not  actually  orthogonal,  whenever  an  exchange  integral  appears  in  the  formulation  we 
set  it  equal  to  a Heisenberg  exchange  integral  J.  The  lowest  singlet  and  triplet  ener- 
gies will  be  determined  as  a function  of  J as  will  the  energy  difference  of  these  lowest 
states.  For  this  model  there  is  just  one  type  of  exchange  integral  appearing  that  being 
between  ar  orbital  located  on  one  of  the  He*  ions  and  the  s orbital  on  the  central  He 
atoiu.  Again  we  restrict  cur  considerations  to  the  superexchange  configurations.  The 
results  for  this  approach  are  given  below. 


Singlets 

T riplets 

QQ  = - 15.  639503  - 2J 

11 

o 

o 

X 

- 15.639503 

- 2J 

Qj  = + 0.  338732 

^^01 

+ 0. 338732 

J J = - 15.  390183 

»11 

- i5. 391925 

- 2J 

If  we  solve  this  problem  and  find  the  final  singlet-triplet  energy  separation  in  terms  of 
J,  and  substiiute  the  value  of  the  Heisenberg  exchange  integral,  which  is 

■^Heis.  = {“l  * “z  * “iz} 

= - 0.  024289,  i = 1 or  4. 
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we  find  Ihe  singlet-triplet  energy  separation  to  be 


E(5  = 0)  - E‘S  = 1)  = - 0.  0 16065 


This  energy  difference  is  to  be  compared  with  the  value  obtained  by  the  proper  treat- 
ment of  the  non-orthogonality  problem,  this  value  b.  ing 

E(S  = 0)  - E(S  = 1)  = - 0.010537 


We  see  that  even  for  a problem  in  which  the  overlap  integral  is  quite  small,  which 
means  that  our  orbitals  are  nearly  orthogonal,  an  error  of  app’.  oximately  50  percent 
in  the  final  energy  separation  is  incurred  by  this  kind  of  a treatment. 

Having  found  the  important  configurations  using  non-orthogonal  orbitals,  the 
problem  can  easily  be  reformulated  in  terms  of  orthogonal  orbitals.  First  we  set  up 
an  s-like  function  about  the  central  He  which  is  orthogonal  to  Uj  and  u,| 


s - 6 ^{u  , + u .) 

UjS'  1 4' 


where  s is  the  non-orthogonal  orbital  used  above.  This  gives 


s = s'  + 6 _( 


UjS'^1  " 


U4). 


We  can  now  replace  the  non-orthogonal  s in  the  states  ij/p  and  relating  to  the  super- 
exchange  configurations  by  its  expansion  in  orthogonal  orbitals.  When  this  is  done  and 
the  determinants  expanded,  it  is  easily  recognized  that  by  using  the  configurations 
given  in  Table  6-4  and  carrying  out  the  configuration  interaction  with  these  orthogonal 
orbitals,  one  would  get  a final  state  of  energy  at  least  as  low  as  that  found  by  using 
non-orthogonal  orbitals. 


Table  6-4 
u,  s' 

1 2 1 

2 1 1 

1 1 2 

2 0 2 

The  same  sort  of  procedure  can  be  used  to  find  the  correct  set  of  configura- 
tions for  the  — ivin^^  system  previously  considered.  For  example,  substituting 

dir  for  Uj,  pa  for  s',  and  d'a  for  u^  in  Table  6-4  shows  the  configurations  which  must 
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be  included  in  order  to  describe  the  transfer  of  a pa  electron  on  the  0~  to  the  da  of  one 
of  the  Mn^^  ions  in  terms  of  orthogonal  orbitals. 

Suppose  we  consider  a treatment  of  the  present  model  by  orthogonal  orbitals, 
but  in  which  wc  IcaVc:  out  the  final  v_onfiguration  given  in  Table  6-4  in  which  both  elec- 

j. 

trons  on  the  central  He  have  been  transferred,  one  to  each  He  ion.  It  can  be  shown 
that  if  this  extra  configuration  is  not  included,  it  is  impossible  to  get  an  antiferromag- 
netic coupling.  To  do  this,  we  proceed  as  in  the  last  Progress  Report  by  setting  up 
the  Hamiltonian  matrix  for  the  singlet  states  and  equating  this  to  the  Hamiltonian  ma- 
trix for  the  triolet  states  olus  a difference  matrix  D. 


»S=0  = "s=l  ^ ° 


The  form  of  D for  the  superexchange  configuration  is 


o 

o 

o 

D 

n 


where  0^2  ts 

D^2  = + ^(UjS|su^)  + 2{u,s|sUj)  {6 

Since  is  positive  definite  the  D matrix  cannot  have  a negative  eigenvalue.  If 
is  the  eigenvector  corresponding  to  the  lowest  lying  energy  Eq  of  Hg  _ g,  then 

^*0^^S  = 0^*0^  “ 
or 

Eq  = <fo,  " D (6, 

But  is  surely  5 Ej  the  lowest  triplet  encrg>-,  hence 

Eg  » Ej  + D (6, 

Knowing  that  the  matrix  D has  no  negative  eigenvalue  shows  us  that 

(♦g|D|*g)  ^ 0 (6, 


i 
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and 

Eq  i Ej  (6.  10) 

Therefore,  the  Icv.cot  singlet  energy  cannot  lie  below  the  lowest  triplet  energy  if  one 
exeludes  the  extra  eonfiguration.  It  is  apparent  that  when  using  orthogonal  orbitals, 
one  must  not  leave  out  any  configuration  that  will  couple  strongly  with  the  low  lying 
states. 


Reference 

1.  r*.  \V.  Aiideison,  Phys.  Rev.  79.  350  (i960). 


G.  W.  Pratt,  Jr. 
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7.  A SYSTEM  OF  TWO  p ELECTRONS 

As  was  described  in  the  previous  Progress  Report,  a system  of  two  p 

electrons  around  two  atomic  cores  is  treated  using  Gaussian  wave  functions  and  taking 

into  account  the  configuration  interactions.  Among  all  the  possible  configurations, 

those  treated  are  show.n  in  Table  7-1.  The  numbers  in  the  table  mean  the  number  of 

i 2,\ 

electrons  in  the  corresponding  molecular  orbiial.  ' ' 

'T'-blc  7-1 


States 

1 

2 

3 

4 

5 

6 


+ + 
(TV  IT 

U g U 


Types  of  configurations  treated  were  and  'tT  , each  of  which  having  con- 

& ^g  g 

tribution  from  states  1,  2,  3 and  4;  3 and  4;  5 and  6,  respectively. 

The  Hamiltonian  has  the  same  form  as  the  one  for  the  hydrogen  molecule. 
The  atomic  wave  function  has  a form,  for  instance 


P 


o 


= 2 yb  (-^y  (z  + exp  - b(iT  + 


z-direction  being  taken  along  the  line  joining  the  two  nuclei,  and  R being  the  internu- 
clear  distance.  The  only  parameter  we  have  is  b,  which  was  chosen  as  unity  in  the 
actual  numerical  computation  in  order  to  obtain  qualitative  information. 

The  secular  equations  for  and  ^TT  are  quadratic,  but  those  for  ^ 

states  are  four  by  four.  It  may  be  of  some  interest  to  compare  the  diagonal  elements 
and  the  eigenvalues  of  four  by  four  matrices.  These  are  shown  in  Fig.  7-1,  which 
looks  fairly  complicated.  Atomic  units  are  used  for  the  internuclear  distances  and  the 
energies. 


In  Fig.  7-2,  the  lowest  energy  values  of  the  three  types  of  states  are  shown. 
From  this  figure,  one  can  conclude  the  following:  (l)  The  lowest  state  is  sometimes 
singlet  and  sometimes  triplet.  (2)  It  is  difficult  to  derive  any  definite  relation  between 
the  present  model  of  two  p-electron  systems  and  ferromagnetism  which  might  be  ex- 
pected for  the  crystal  composed  of  p-electron  atoms. 

The  author  wishes  to  thank  Gloria  Freedman  for  the  numerical  computation 
of  the  four  by  four  secular  equations. 
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Fig.  7-2 


(A  SYSTEM  OF  TWO  p ELECT KONS) 
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H.  Kikuchi 


8.  CONFIGURATIONAL  INTERACTION  APPLIED  TO  THE  HYDROGEN  MOLECULE 


In  the  molecular  orbital  treatment  of  molecules  one  commonly  constructs  the 
orbitals  out  of  orthogonalized  linear  combi  nations  of  atomic  orbitals.  By  choosing  a 
sufficiently  large  number  of  these  functions,  one  can  get  an  arbitrarily  good  representa- 
tion of  the  true  wave  function.  This  scheme  has  one  unhappy  aspect,  however.  As  the 
difference  between  the  total  number  of  orbitals  employed  and  the  number  occupied  in 
the  isolated  atoms  is  increased,  the  number  of  configurations  grows  at  a prodigious 
rate,  leading  to  secular  equations  of  very  great  size.  One  may  hope  that  by  employing 
molecular  orbitals  constructed  more  accurately  than  as  linear  combinations  of  atomic 
orbitals  one  can  get  a good  representation  of  the  many -electron  wave  function  with  only 
a small  number  of  configurations. 

With  this  hope  in  mind  various  molecular  orbitals  are  being  used  in  a study  of 
the  hydrogen  molecule.  One  can  approximate  the  true  Hamiltonian  by  the  spin  free 
terms,  kinetic  energy  plus  electrostatic  interactions  between  the  four  particles 


H = - V,"  - v/ 


2 ^ 2* 
r, , R 


In  this  case  one  may  factor  out  the  antisymmetric  spin  factor  and  look  for  a solution  of 
the  form 

‘^'(1,2)  = ^ I ^ {aU)P(2)  - P(l)a(2)j. 

i Jk 

where  the  o-'s  are  one-electron  molecular  orbitals.  Each  term  in  these  summations 
is  a singlet,  and  must  be  so  chosen  as  to  satisfy  the  symmetry,  ^ of  the  observed 
electronic  structure. 

As  the  first  configuration  I chose  a function  (l)  (2)  constructed  from 

( I ) ^ ° 

molecular  orbitals  determined  by  Coulson 

o-  =Ne  (1  + hji  + k\) 


where  >.  and  p are  the  elliptic  coordinates 


‘a  " ‘B 


‘ a ’’b 


For  the  internuclear  separation  R = 1.4  a.  u.  Coulson  found  h and  k by  minimizing 


The  13.  54  e.  v.  atomic  units  are  used  here. 
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<o-^_(l)  (t’  (2)lHla''^  (1)  o-^  (2)> 

E(h.k)  = i- g-  i S 

<'  s'"  B> 

It  seems  desirable  to  add  to  (r^p(l)  (r^g{2)  another  configuration  made  of  one- 
electron  functions  antisymmetric  in  the  nuclei.  Coulson's  symmetric  orbiials  are  ap- 
proximate solutions  of  the  self-consistent  field  equation 

— + V )4>  = £4> 

^A  Tg  H 

where 

V (1)  = fff*  ’^(2)-^/  (2)  dr. 

As  an  antisymmetric  orbital,  an  approximation  to  the  lowest  antisymmetric  solution  of 
this  same  Sohrddinger  equation  has  been  constructed  in  the  form 

= N p e‘  + cp^  + d\) 

adjusting  c and  d to  minim' ze  the  energy.  Having  thus  found  and  , a configura- 

& ^ + 

tion  interaction  was  set  ud  between  the  configurations  when  both  electrons  were  in  a , 

+ “ 
and  where  both  were  in  <r  The  resulting  secular  equation,  which  has  been  solved,  is 


+ I + 2 

^(7  (7  > 

^ U ' U^ 


The  integrals  involved  are  all  trivial  except  for  the  j — integrals,  which  can  be  very 
messy.  But  as  one  is  able  to  approximate  solutions  in  the  form  of  exponentials  n X. 
and  ^ times  a polynomial  in  vof-inhios.  ♦*’«•  Rudenberg^^^  scheme  may  be  used 

directly,  and  proves  very  convenient.  The  papers  of  Kotani  et.  al.  give  a partial 
listing  of  a set  of  functions  needed  in  the  integration.  This  table  is  presently  being 
extended  by  C.  C.  J.  Roothaan  and  Klaus  Riidenberg.  In  the  above  quadratic  secular 
equation  the  elements  are 
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(2)lH|  ff  (l)  o-^  (2)> 

2 U B S 

^ g' 


= - 2.  26 1 37  a.  u. 


(l)  (2)1H|  o-^, ,{2)> 


u u 


u u 


1/  y 


. 16778  a.  u. 


<]ff^  lo-^  > ^ 

g g^  ^ u \ir 


. 26528  a.  u. 


Binding  Energy 
- 3.  535  e.  V. 


29.  48  e.  V. 


The  binding  energy  resulting  from  the  configuration  interaction  is  - 3.944  e.  v.  , com- 
pared to  the  Coulson  values  - 3.  535  e.  v. , an  increase  of  - . 409  e.  v.  by  the  addition 
of  the  second  configuration,  which  is  itself  energetically  very  high.  The  experi- 
mentally observed  binding  energy  is  - 4.  72  e.  v. 

There  are  still  a number  of  interesting  calculations  to  be  made  in  this  prob- 
lem. If  the  energies  of  the  two  configurations  were  closer  together,  the  lower  would 

be  pushed  down  more.  Given  <r  , and  the  form  of  a , the  values  of  the  coefficients  in 

g 

<r^  which  minimize  the  total  energy  are  the  best  ones  to  use,  and  these  may  not  be  ex- 
actly the  same  as  the  values  which  minimize  the  energy  of  the  average  field  equation. 

I am  now  trying  to  solve  the  secular  equation  with  the  parameters  undetermined,  to 
see  how  well  one  can  do  with  this  configuration,  and  how  far  this  solution  is  from  the 
solution  of  the  Hartree  equation. 
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9.  THE  WATER  MOLECULE 

Since  the  last  Quarterly  Progress  Report  we  have  completed  all  the  one-center 
electrostatic  integrals  involved  in  the  calculation  of  the  energy  for  the  equilibrium  in 
ternuclear  distance  of  1.  8 a.u.  We  have  also  completed  many  of  the  two-center  inte- 
grals including  the  exchange  integrals  involving  the  hydrogen  Is  and  the  oxygen  Is  and 
2s  wave  functions.  These  expressions  involve  infinite  sums  arising  from  the  expan- 
sion of  the  hydrogen  wave  function  about  the  oxygen  center.  We  have  carried  six  terms 

in  the  Legendre  function  expansion  (through  Pc),  with  the  result  that  the  neglected  Is 

-6  ^ 

terms  contribute  only  about  10  of  the  total  integral,  while  the  error  in  the  2s  sum 
(after  a simple  extrapolation)  was  of  the  order  of  10 

G.  F.  Koster  and  H.  C.  Schweinler 


10.  MECHANIZA'IION  OK  MOLECULAK  ORBITAL  DETERMINATION 


The  Whirlwind  program  for  the  self-consistent  determination  of  the  optimum 
linear  combinations  of  atomic  or’ntals  is  in  the  testing  and  assembling  stage.  The 
complete  p:  ocedure  has  been  broken  down  into  smaller  entities,  much  like  the  flow 
chart  shown  in  the  last  Progress  Report,  and  each  part  is  being  made  error-free  pre- 
paratoj  y to  combining  irto  cne  process  continuous  from  the  input  of  integrals  to  the 
output  of  self-consistent  eigenvectors  and  energy  levels. 

A.  Meckler 


11.  LIMITED  CONFIGUKAIION  INTERACTiOW  1 KKATMENT  OF  THE  NH^  MOLECULE 

The  greatest  portion  of  the  analytic  work  described  in  the  last  Progress  Re- 
port has  been  completed.  Numerical  calcalation  of  all  necessary  integrals  has  been 
started  for  the  observed  internuclear  separations. 

In  order  tc  test  the  accuracy  of  the  numerical  integration  scheme,  which  uses 
a fine  mesh  to  compensate  for  the  relative  crudity  of  Simpson's  Rule,  the  overlap  in- 
tegral between  two  hydrogen  Is  functions  on  different  centers  was  computed  for  the 
internuclear  separation  of  3 a.  u.  and  the  result  compared  with  the  exact  answer.  The 
results  are  given  below. 


Numerical  Integration  0.  3485156 
Exact  Answer  0.  3485095 

Difference  0.0000061 

Although  the  hydrogen  overlap  integral  is  not  the  most  stringent  test  of  the  method  of 
integration,  the  very  accurate  result  obtained  plus  the  fact  tnat  the  final  energy  should 
be  in  the  first  decimal  place  seem  to  warrant  the  conclusion  that  the  accuracy  of  the 
final  result  will  not  be  limited  by  the  scheme  of  integration. 

II.  Kaplan 
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!2.  CON FIGU RATION  INTERACTION  FOR  THE  FLUORINE  MOLECULE 


(1) 


By  choosing  analytic  wave  functions,  it  was  hoped  that  the  many  formulas  and 
tables  developed  for  the  purpose  could  be  used  to  evaluate  the  two-center  atomic  in- 
tegrals necessary  for  the  problem.  Kopineck's^^^  tables,  based  on  the  work  of  Kotani, 
Amemiya,  and  Simose,  are  in  the  most  convenient  form  for  use.  However,  these 
tables  contain  only  the  integrals  involving  just  the  2s  and  2p  functions  on  the  two  a*oms. 
The  remainder  of  the  integrals  can  be  computed  by  means  of  formulas.  Part  of  this 
remainder  can  be  expressed  in  terms  of  simple  expressions,  but  many  can  only  be  ex- 
pressed in  very  complicated  forms.  Because  of  the  limited  time  that  the  author  had 
with  the  Group,  it  was  not  possible  to  complete  the  necessary  calculations.  Indeed  the 
time  was  so  short  that  it  did  not  even  seem  expedient  to  start  the  numerical  computa- 
tion. At  the  present  time  the  author  cannot  say  anything  about  the  future  progress  of 
the  problem. 
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13.  NUCLEAR  ELECTRIC  QUADRUPLE  INTERACTION  IN  THE  KCl  MOLECULE 


Investigation  of  the  polarized  one-electron  function; 

4«  - + Cr‘'  cos  e) 

where  = unperturbed  atomic  function 

C,  V = variation  parameters 

is  continuing.  A calculation  using  these  one-electron  functions  in  a single  determinant 
representation  of  F in  a uniform  electric  field  is  in  progress. 

L.  C.  Allen 


14.  SCATTERING  OF  NEUTRONS  BY  O^^ 


Calculation  of  the  slow  neutron  differential  scattering  cross-section  for  O, 

( 1 ) ^ 

IS  continuing  along  the  lines  indicated  in  a previous  Report.  ' The  cross-section 
takes  account  of  inelastic  transitions  between  rotational  levels  of  the  molecule.  Effort 
IS  now  being  concentrated  on  the  evaluation  of  a number  of  integrals  which  arise  in 
connection  with  the  scattering  associated  with  the  magnetic  spin-spin  part  of  the  neutron- 
electron  interaction.  It  is  planned  to  obtain  a total  cross-section  for  comparison  with 
results  of  experiments  under  way  at  Brookhaven  National  Laboratory. 
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15.  TRANSITION  PROBABILITIES  FOR  X-RAY  EMISSION  IN  THE  TRANSITION  METALS 


One  of  the  most  direct  methods  used  to  obtain  the  characteristics  of  the  en- 
ergy distribution  of  the  conduction  electrons  in  a metal  is  through  the  study  of  their 
soft  x-ray  emission  spectra.  This  radiation  is  caused  by  the  transition  of  an  elec- 
tron from  the  continuous  band  of  filled  valence  levels  into  an  empty  inner  level  of  an 
atom  ionized  by  electron  impact.  The  interaction  between  the  lattice  potential  and 
these  inner  atomic  levels  is  weak,  and  therefore  the  corresponding  energy  levels  re- 
main fairly  sharp.  The  continuous  band  of  radiation  will  therefore  be  a measure  of 
the  distribution  of  the  occupied  valence  electrons.  The  correlation  cannot  however  be 
made  directly  between  the  intensity  of  the  emitted  radiation  and  the  density  of  occu- 
pied energy  states  because  the  observed  intensity  depends  on  the  transition  probability 
between  states.  If  dN/dE  denotes  the  number  of  electrons  per  unit  energy  range  and 
f(E)  the  probability  that  an  electron  in  an  energy  state  E will  make  a transition  to  an 
inner  level,  the  intensity,  l(E),  is  given^*’^^  as 

1(E)  = f(E)  dN/dE  (15.  1) 

In  the  transition  metals  the  valence  band  is  derived  from  the  4s,  4p,  3d  and 

higher  levels,  and  is  pictured  in  band  theory  as  a wide  low-density  s-p  band  capable 

of  holding  2 electrons  per  atom  with  a rather  narrow  d band,  capable  of  holding  10 

electrons  per  atom,  overlapping  it.  Because  of  the  higher  density  of  states  in  the 

d band  the  observed  emission  bands  would  be  expected  to  represent  only  the  five  d 
(4) 

bands.  Gyorgy  has  recently  measured  the  emission  spectra  from  the  transition 
metals  and  has  found  a considerable  contribution  from  the  s-like  states  with  a cor- 
respondingly lower  intensity  from  the  d bands.  In  order  to  explain  these  and  other 
results  he  is  led  to  assume  that  the  transitions  from  the  s-p  band  to  the  3P  level  have 
a greater  probability  than  transitions  from  the  d bands  to  the  same  atomic  level.  The 
probability  of  a transition  of  an  electron  from  a valence  state  + to  an  inner  level  4> 
is  proportional  to 


P = ^ j ^ + ^y<})*z4«iTj  ' 


(15.2) 


.•Hid  depends  on  the  exact  nature  of  the  wave  functions  involved.  Since  these  wave 
functions  are  not  known  exactly,  any  calculation  must  therefore  give  p only  approxi- 
mately. For  the  large  effect  we  are  seeking,  even  an  order  of  magnitude  calculation 
seems  to  be  instructive  at  this  time. 

For  the  3d  ;ind  more  inner  levels  we  choose  as  our  crystal  wave  functions  a 
spatial  Bloch  sum  of  atomic  wave  functions,  ^(r^ 


>l<j.(r)  = Z e " - r ) 

K n " 


(15.  3) 
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where  <l)(r  - r ) is  the  atomic  wave  function  centered  at  r , the  vector  from  the  origin 

n'  th  ” 

of  coordinates  to  the  center  of  the  n unit  cell.  The  4s  wave  function  is  chosen  as  a 

(5) 

plane  wave  orthogonalized  ' to  tlic  inner  Bloch  wave  functions  defined  above.  The 
quantity  p is  then  computed  fro».i  Eq.  (15.2)  for  the  case  when  the  initial  state  is  3d 
and  compared  with  the  value  obtained  for  an  initial  4s  state.  We  have  fitted  the  Har- 
tree^^^  one-electron  wave  functions  for  chromium  by  means  of  a polynoniial  times  a 
sum  of  exponentials  in  order  to  obtain  analytic  expressions  for  our  <j>(r).  Chromium 
was  chosen  because  it  is  the  element  for  which  Gyorgy  has  found  the  emission  bands 
and  for  which  self-consistent  field  wave  functions  arc  available. 
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